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We present a new approach to describe tlie dynamics of an isolated, gravitationally bound as- 
tronomical A'^-body system in the weak field and slow-motion approximation of general relativity. 
Celestial bodies are described using an arbitrary energy-momentum tensor and assumed to possess 
any number of internal multipole moments. The solution of the gravitational field equations in any 
reference frame is presented as a sum of three terms: i) the inertial flat spacetime in that frame, 
ii) unperturbed solutions for each body in the system boosted to the coordinates of this frame, and 
iii) the gravitational interaction term. Such an ansatz allows us to reconstruct all features of the 
gravitational field and to develop a theory of relativistic reference frames. We use the harmonic 
gauge conditions to impose a significant constraint on the structure of the post-Galilean coordinate 
transformation functions that relate global coordinates in the inertial reference frame to the local 
coordinates of the non-inertial frame associated with a particular body. The remaining parts of 
these functions are constrained using dynamical conditions, which are obtained by constructing the 
relativistic proper reference frame associated with a particular body. In this frame, the effect of 
external forces acting on the body is balanced by the fictitious frame-reaction force that is needed 
to keep the body at rest with respect to the frame, conserving its relativistic linear momentum. 
We find that this is sufficient to determine explicitly all the terms of the coordinate transformation. 
The same method is then used to develop the inverse transformations. The resulting post-Galilean 
coordinate transformations have an approximate group structure that extends the Poincare group 
of global transformations to the case of a gravitational A^-body system. We present and discuss the 
structure of the metric tensors corresponding to the reference frames involved, the rules for trans- 
forming relativistic gravitational potentials, the coordinate transformations between frames and the 
resulting relativistic equations of motion. 

PACS numbers: 03.30.-|-p, 04.25. Nx, 04.80.-y, 06.30.Gv, 95. 10. Eg, 95.10.Jk, 95.55.Pe 

I. INTRODUCTION 

Recent experiments have successfully tested Einstein's general theory of relativity in a variety of ways to remarkable 
precision [l|, |2|. Diverse experimental techniques were used to test relativistic gravity in the solar system, namely: 
spacecraft Doppler tracking, planetary ranging, lunar laser ranging, dedicated gravity experiments in space and many 
ground-based efforts [2tJ]- Given this phenomenological success, general relativity became the standard theory of 
gravitation, especially where the needs of astronomy, astrophysics, cosmology and fundamental physics are concerned 
y. The theory is used for many practical purposes involving spacecraft navigation, geodesy and time transfer. It is 
used to determine the orbits of planets and spacecraft and to describe the propagation of electromagnetic waves in 
spacetime [2[. 

As we shall see, finding a solution to the Einstein's equations in the case of an unperturbed one body problem is 
quite a simple task. However, it turns out that a generalization of the resulting post-Newtonian solution to a system 
of N extended arbitrary bodies is not straightforward. 

A neutral point test particle with no angular momentum follows a geodesic that is completely defined by the 
external gravitational field. However, the coupling of the intrinsic multipole moments of an extended body to the 
background gravitational field (present due to external gravitational sources), affects the equations of motion of such 
a body. Similarly, if a test particle is spinning, its equations of motion must account for the coupling of the body's 
angular momentum to the external gravitational field. As a result, one must be able to describe the interaction of 
a body's intrinsic multipole moments and angular momentum with the surrounding gravitational field. Multipole 
moments are well-defined in the local quasi-inertial reference frame generalizing what was defined for the unperturbed 
one-body problem. While transforming these quantities from one coordinate frame to another, one should account 
for the fact that the gravitational interaction is non-linear and, therefore, these moments interact with gravitational 
fields, affecting the body's motion. 

When the Riemannian geometry of the general theory of relativity is concerned, it is well known that coordinate 
charts are merely labels. Usually, spacetime coordinates have no direct physical meaning and it is essential to 



construct the obscrvablcs as coordinate-independent quantities. Thus, in order to interpret the results of observations 
or experiments, one picks a specific coordinate system, chosen for the sake of convenience and calculational expediency, 
formulates a coordinate picture of the measurement procedure, and then one derives the observable out of it. It is 
also known that an ill-defined reference frame may lead to appearance of non-physical terms that may significantly 
complicate the interpretation of the data collected |5| . Therefore, in practical problems involving relativistic reference 
frames, choosing the right coordinate system with clearly understood properties is of paramount importance, even as 
one recognizes that in principle, all (non-degenerate) coordinate systems are created equal [6|. 

Before one can solve the global problem (the motion of the entire A^-body system), the local gravitational problem 
(in the body's vicinity) must be solved first. The correspondence of the global and local problems is established using 
coordinate transformations by which representations ofphysical quantities are transformed from the coordinates of 
one reference frame to another. It is understood (see [7J, |8| for discussion) that, in order to present a complete solution 
for TV-body problem in the general theory of relativity, one must therefore find the solutions to the following three 
intertwined problems: 

1). The global problem: We must construct a global inertial frame, such as a barycentric inertial reference frame, for 
the system under study. This frame can be used to describe the global translational motion of N extended bodies 
comprising the system; 

2). The local problem: For each body in the system, we must construct a local reference frame. This frame is used to 
study the gravitational field in the vicinity of a particular body, to establish its internal structure and to determine 
the features of its rotational motion; and finally, 

3). A theory of coordinate reference frames: We must establish rules of coordinate transformations between the global 
and local frames and use these rules to describe physical processes in either of the two classes of reference frames. 

A prerequisite to solving the first two of these problems is knowing the transformation rules between global and 
local reference frames. Thus, establishing the theory of astronomical reference frames and solving the equations of 
motion for celestial bodies become inseparable. 

Modern theories of relativistic reference frames 0, IMtIH , dealing predominantly with the general theory of relativity, 
usually take the following approach: As a rule, before solving gravitational field equations, four restrictions (coordinate 
or gauge conditions) are imposed on the components of the Riemannian metric gmn- These conditions extract a 
particular subset from an infinite set of spacctime coordinates. Within this subset, the coordinates are linked by 
smooth differentiable transformations that do not change the coordinate conditions that were imposed. A set of 
differential coordinate conditions used in leading theories of relativistic reference systems, such as that recommended 
by the International Astronomical Union (see, for instance [6|), are the harmonic gauge conditions. In addition, a set 
of specific conditions designed to fix a particular reference frame is added to eliminate most of the remaining degrees 
of freedom, yielding an explicit form for the coordinate system associated with either frame. 

In a recent paper [l6|, we studied accelerating relativistic reference frames in Minkowski spacetime under the 
harmonic gauge condition. We showed that the harmonic gauge, which has a very prominent role in gravitational 
physics starting with the work of Fock [l3| , allows us to present the accelerated metric in an elegant form that depends 
only on two harmonic potentials. It also allows reconstruction of a significant part of the structure of the post-Galilean 
coordinate transformation functions relating inertial and accelerating frames. In fact, using the harmonic gauge, we 
develop the structure of both the direct and inverse coordinate transformations between inertial and accelerated 
reference frames. Such a complete set of transformations cannot be found in the current literature, since usually 
either the direct [a, |ll| or the inverse [lO| transformation is developed, but not both at the same time. The remaining 
parts of these functions were determined from dynamical conditions, obtained by constructing the relativistic proper 
reference frame of an accelerated test particle. In this frame, the effect of external forces acting on the observer are 
balanced by the fictitious frame-reaction force that is needed to keep the test particle at rest with respect to the frame, 
conserving its relativistic linear momentum. We find that this is sufficient to determine explicitly all the terms of 
the coordinate transformation. The same method is then used to develop the inverse transformations. The resulting 
post-Galilean coordinate transformations exhibit an approximate group structure that extends the Poincare group of 
global spacetime transformations to the case of arbitrarily accelerated observers moving in the gravitationbal field. 

In the present paper we continue the development of this method and study the dynamics of a gravitationally 
bounded astronomical A^-body system. To constrain the available degrees of freedom, we will again use the harmonic 
gauge. Using the harmonic gauge condition, we develop the structure of both the direct and inverse harmonic 
coordinate transformations between global and local reference frames. The method presented in this paper could 
help to generalize contemporary theories of relativistic reference frames and to allow one to deal naturally with both 
transformations and the relevant equations of motion in a unified formalism. Finally, the method we present does 
not rely on a particular theory of gravitation; instead, it uses covariant coordinate transformations to explore the 



dynamics in the Minkowski spacetimc from a general perspective. Thus, any description from a standpoint of a metric 
theory of gravity must have our results in the general relativistic limit. 

The outline of this paper is as follows: 

In Sec. In] we introduce notation and briefly study the trivial case of a single unperturbed body, followed by a 
description of the motion of a system of N weakly interacting, self-gravitating, deformable bodies by generalizing 
the approximate Lorentz transformations. This is accomplished by introducing acceleration-dependent terms in the 
coordinate transformations. The resulting transformation is given in a general form that relies on a set of functions 
that arc precisely determined in the subsequent sections. 

The local coordinate system of an accelerated observer is not unique. We use the harmonic gauge to constrain 
the set of co-moving coordinate systems in the accelerated reference frame. We carry out this task and determine 
the metric tensor describing the accelerated reference frame and the structure of the coordinate transformations that 
satisfy the harmonic gauge conditions. We observe that the metric in the accelerated frame has an elegant form that 
depends only on two harmonic potentials, which yield a powerful tool that allows the reconstruction, presented in 
Sec. mil of a significant part of the structure of the post-Galilean coordinate transformation functions between inertial 
and accelerating frames. 

To fix the remaining degrees of freedom and to specify the proper reference frame of an accelerated observer, 
we introduce a set of dynamical conditions in Section IIVI Specifically, we require that the relativistic linear three- 
momentum of the accelerated observer in its proper reference frame must be conserved. This conservation introduces 
the necessary additional constraints to determine uniquely all remaining terms. 

A similar approach can also be carried out in reverse, establishing coordinate transformation rules from an accel- 
erating to an inertial reference frame, which is done in Sec. |Vl Key to the approach presented in this section is the 
use of the contravariant metric for the accelerating frame, which leads in a straightforward manner to the inverse 
Jacobian matrix, allowing us to present the inverse transformations in which the roles of the inertial and accelerating 
coordinates are reversed. Although the calculations are formally very similar to those presented in Sees. Ill Dl and HVl 
subtle differences exist; for this reason, and in order to keep Sec.|V]as self-contained as possible, we opted in favor of 
some repetitiveness. 

We conclude by discussing these results and presenting our recommendations for future research in Section FVIl We 
also show the correspondence of our results to those obtained previously by other authors. 

II. MOTION OF AN ISOLATED AT-BODY SYSTEM 

Our first goal is to investigate the case of N bodies, forming a gravitationally bound system that is isolated, i.e., 
free of external gravitational fields, described by a metric that is asymptotically fiat. Thus, sufficiently far from the 
iV-body system the structure of its gravitational field should resemble that of the one-body system. 

A. The unperturbed one-body system 

From a theoretical standpoint, general relativity represents gravitation as a tensor field with universal coupling to 
the particles and fields of the Standard Model. It describes gravity as a universal deformation of the flat spacetime 
Minkowski metric, jmn' 

Alternatively, it can also be defined as the unique, consistent, local theory of a massless spin-2 field hmn, whose source 
is the total, conserved energy-momentum tensor (see [la 13 ^'^'^ references therein).^ 



1 



The notational conventions employed here are those used by Landau and Lifshitz 1201 : Letters from the second half of the Latin alphabet, 
m,n,... = 0...3 denote spacetimc indices. Greek letters a,l3,... = 1...3 denote spatial indices. The metric 7mn is that of Minkowski 
spacetime with 7mn = diag(-|-l, — 1, — 1, — 1) in the Cartesian representation. The coordinates are formed such that {c4, r) = (x^'jX"), 
where c is the velocity of light. The semicolon a;rn = Vm denotes a covariant derivative. We employ the Einstein summation convention 
with indices being lowered or raised using fmn ■ Round brackets surrounding indices denote symmetrization and square brackets denote 
anti-symmetrization. A dot over any function means a differentiation with respect to time t, defined by x^ = ct. We use negative powers 
of c as a bookkeeping device for order terms; g^lj, fc = 1, 2, 3... denotes the fc-th term in the series expansion of the metric tensor. Other 
notations will be explained in the paper. 



Classically [2l|, |22| , the general theory of relativity can be defined by postulating the action describing the gravita- 
tional field and its coupling to matter fields. Absent a cosniological constant, the propagation and self-interaction of 
the gravitational field is described by the action 



Scigmn] = 777-7=; / d^Xy/^R, (2) 




where G is Newton's constant, g™" is the inverse of gmm 9 = det^mn, R is the trace of the Ricci tensor that in turn 
is given by R^n = dhT'^^ - ^m^nk + ^mn^ii ~ ^nii^nk' and rf„„ = \g^'P{dm.gpn + dnQprn - dpSmn) are the Christoffel 
symbols. 

The coupling of g,„„ to all matter fields (this would generally mean all the fields of the Standard Model of particle 
physics) is accomplished by using it to replace the Minkowski metric everywhere [3|. Varying the total action 

Stotli^, AmiH] gmn\ ~ SQ[gmn] + '5sm[V'j ^m, H ] gmn], (3) 

with respect to gmn we obtain Einstein's field equations for gravity: 

R =^[T --g T), (4) 

where T"*" = {2/,/^)SCsM/Sgmn is the (presumed to by symmetric) energy-momentum tensor of matter, represented 
by the Lagrangian density £sm- As the density of the Einstein's tensor is conserved, namely S/kW—g{R"^'' — ^g"^''R)] = 
0, it follows that the density of the energy-momentum tensor of matter y^—gT'^^ is also conserved, obeying the 
following covariant conservation equation: 



^mk\ 



0. (5) 



Einstein's equations @ connect the geometry of a four-dimensional Riemannian manifold representing spacetime to 
the stress-energy-momentum of matter contained in that spacetime. The theory is invariant under arbitrary coordinate 
transformations: a;'™ = /"(a;"). This freedom to choose coordinates allows us to introduce gauge conditions that 
may help with solving the field equations ^. For instance, in analogy with the Lorenz gauge (9m ^"^ = 0) of 
electromagnetism, the harmonic gauge corresponds to imposing the condition [l7|: 

5„(y^g"")=0. (6) 

We can now proceed with finding a solution to the Einstein's equations (|4|) that satisfy the harmonic gauge ([6]). 

B. Solution to the unperturbed one-body problem 

Equations (J4|) represent a set of non-linear hyperbolic-type differential equations of the second order with respect to 
the metric tensor of the Riemannian spacetime. This non-linearity makes finding a solution to this set of equations in 
the general case to be a complicated problem to which no analytical solution is known and a full numerical treatment is 
needed. Depending on a particular situation one usually introduces relevant small parameters and develops a solution 
iteratively. 

When studying a problem in the weak gravitational field and slow motion approximation, one often uses the ratio 
v/c of the velocity of motion of the bodies v with respect to baryccntric coordinate reference frame to the speed of 
light c. For the bodies of the solar system this ratio is a small parameter and typically is of the order v/c « 10~^. 
In our formulation, we use the dimensioned parameter c~^ as a bookkeeping device for these order terms. Since we 
are considering a gravitationally bound TV-body system in the weak field and slow motion approximation, there exist 
relations linked by the virial theorem: v"^ /c^ ~ GM/{c^R) [Ij, 1231, with M, R, and v being the mass, radius and 
escape velocity of the body. Thus, when we write j/o + c~^/C, for instance, this implies that K, is of order w^y° with 
V <^ c being the velocity of an accelerating observer or test particle under consideration and c~^/C is of order (w/c)^y°, 
which remains small relative to y". 

To find the solution of the gravitational field equations in first post-Newtonian approximation of the general theory 
of relativity, we expand the metric tensor gmn (HJ '■ 

500 = 1 + c-'g^oo + c-'g^oo + 0{c-'), (7) 

goc = C-3g[3j+0(c-5), (8) 



gcp = icp + c-^gi^l + 0(c-4), (9) 

where 7q^ = (— 1,— 1,— 1) is the spatial part of the Minkowski metric 7mn- Note that goa starts with a term of order 
c~^ , which is the lowest order gravitational contribution to this component. 

In our calculations we constrain the remaining coordinate freedom in the field equations by imposing the covariant 
harmonic gauge conditions in the following form of Eq. ([6]). Thus, for n = and n = a, we correspondingly have: 

Id^il'-gf} - 5|,o') - d'^il = 0(c-5), (10) 

\d\i^+l'''9f})~l'''d'^9^^l = 0{c-\ (11) 

The metric tensor given by Eqs. ([T])-® and the gauge conditions given by Eqs. (fT0 |) -([TT 1) allow one to simplify the 
expressions for the Ricci tensor and to present its contravariant components /J'"" — g"''^ g^''^ Rki in the following form: 

i?"" = c-3iAg[3]0"+o(c-5), (13) 

i?"^ = c-2iA(7[2l"/3 + 0(c-4), (14) 

where D = 7™"(9,„5„ = d^Q+d^d^ and A = d^d'^ are the d'Alembert and Laplace operators of the Minkowski spacetime 
correspondingly. 

To solve the gravitational field equations (|3]), we need to specify the form of the energy-momentum tensor. However, 
it is sufficient to make only some very general assumptions on the form of this tensor that are valid within the post- 
Newtonian approximation. Indeed, we will only assume that the components of the energy- momentum tensor, T™", 
have the following form; 

T™ = c2('TMo" + c-2r[2]oo + 0(c-4)\ r""=c(TW"" + 0(c-2)y T"'3=rPl"'3 + 0(c-2). (15) 

This form is sufficient to construct an iterative solution scheme in the post-Newtonian approximation of the general 
theory of relativity leaving the precise form of the energy-momentum tensor T™" unspecified. These assumptions 
allow us to present the "source term" , 5""" = T™" — ^g™"-T, on the right-hand side of Eq. (|4]) in the following form: 

S"" = lc'{a + 0{c-')), (16) 

5"^ = -7"''^c2a(l + c-22g[2])+^"/3+ 0(^-2). (18) 



where we have introduced the scalar, vector, and shear densities a, cr", and a"": 

a = r[0l00 + C-2(T[2l00-7^;,T[2lA'A)+O(c-4)=^(^T00_^^^y^A^^(^(^-4)^ (^g) 

^a ^ y[l]Oa ^ ^-2j,[3]0a _^ ^^(^-4) ^ lj.Qa _^ 0(^-4)^ (20) 

c 

^a/? ^ y[2]a/3_^a/3^^^y[2]pA^(^(^-2) ^ j,a/3 _ ^a^^^^y^A _^ 0(^-2) ^ (2^) 

We can now express the stress-energy tensor T™" in terms of fr™" introduced by Eq. (fT9|) - ((2T|) : 

T"'> ^ c^a-^-f^xa^-^+Oic-^), T^" = c(a" + ^(c-^)) , T"^ = ^"^ ^ 7"^^7MAf^''^ + ^(c"')- (22) 

Substituting the expressions (fT2)) - (rT4l) for the Ricci tensor together with the source term 5™" given by Eqs. (fT6|) - 
pTjl into the field equations (J4]) of the general theory of relativity, we obtain the following equations to determine the 
components of the metric tensor in the post-Newtonian approximation: 

D{^ff'"°" + ^(5'''""-5(5'''°°)')} + ^5,9.g[2l°"(5[2l^^+7^V'l°°) = ^a + Oic~% (23) 



JLa^PIO" = l^^"+0(c-5), (24) 

}-^gm^P ^ _^"/^^^ + 0(c-4). (25) 



We assume that spacetime is asymptotically flat (no external gravitational field far from the body) and there is 
no gravitational radiation coming from outside the body. In terms of perturbations of the the Minkowski metric 
hmn = 5mn — 7mn introduced by Eq. ([1]), the corresponding two boundary conditions |17| have the form 

[{rhran),r + {rhran) fi] = 0. (26) 

Multiplying Eq. (f23| by 7"'^ and summing the result with Eq. (fSS]) allows one to determine (;P]"/3 4--^"^g[2]oo _ 0{c~'^), 
which was derived using the boundary conditions ([26|) . As a result, the system of equations (|23 |) -([25 |) can cquivalently 
be re-written in the following form: 



lim hmn 


= n 


and 


lim 


r^oo 






r— s-oo 


i+r/c— const 






*+r/c— const 



^{^J'''' + ^M''' - Wn"}) = ^-^- + 0{c-% (27) 

4.A9[3]o« - l^a"+0(c-5), (28) 

Solution to this system of equations may be given as below: 

5"° = 1 + 1^ + ^ + 0{c-% r = ^ + 0{c-% 5"^ =: 7"'' - 7-^1^ + 0{c-% (30) 

where the scalar and vector gravitational potentials w and w" are determined from the following harmonic equations: 

Du; = 47rGcr + C'(c-'*), Aw" = 47rGcr" + ©(c^^), (31) 

which, according to the harmonic gauge conditions (|10 p - (|lip and the energy-momentum conservation equation ([5|). 
satisfy the following two Newtonian continuity equations 

cdQW + d,w^ = Oic-"^) and cd^a ^ d^a' ^ OicT"^). (32) 

Assuming that spacetime asymptotically flat (j26p . one can write a solution for w and w°^ in terms of the advanced 
and retarded potentials. The recommended solution Q, half advanced and half retarded, reads 

«;(i,x) = Gy'dV^^ + ^G^y'dVa(t,x')|x-x'|+0(c-3), (33) 

u;"(i,x) = G /dV^^^ 
J |x-x' 



'3-''^"^*''''^+0(c-2). (34) 



Finally, we can present the unperturbed solution for an isolated one-body problem in terms of Minkowski metric 
perturbations /i„j„ introduced by Eq. ([IJ as below: 

9mn — 7mn ~r '^vim 

C - -|^ + 5^ + 0(c-% C = -,J-^ + O(c-), /.i°j = ,J-^ + 0(c-% (35) 



Equations (|35j) represent the well-known solution for the one-body problem in the general theory of relativity 
UIj S UM- The new method, presented here, relies on the properties of this solution in developing a perturbation 
theory needed to flnd a solution in the case of the iV-body problem. 

C. Ansatz for the Af-body problem 

To describe the dynamics of the A^-body problem wc introduce A^+ 1 reference frames, each with its own coordinate 
chart. We need a global coordinate chart {a;'"} defined for an inertial reference frame that covers the entire system 



under consideration. In the immediate vicinity of each body b in the A''-body system, wc can also introduce a chart 
of local coordinates {y™} defined in the frame associated with this body. In the remainder of this paper, we use {a;™} 
to represent the coordinates of the global inertial frame and {j/™} to be the local coordinates of the frame associated 
with a particular body b. 

First, we assume that the studied A'^-body system is isolated and, similarly to the one body problem, there exists 
a barycentric inertial system (with global coordinates denoted as x'' = (ci,x)) that asymptotically resembles the 
properties of the one-body problem discussed above, namely the metric tensor, gmni^'^), describing the A'^-body 
system in that reference frame is asymptotically flat (no gravitational fields far from the system) and there is no 
gravitational radiation coming from outside the system. 

Therefore, in terms of the perturbations of the Minkowski metric ft.™" = g™" — -y™", the boundary conditions are 
identical to those given by (|26p and, in the case of the A^-body problem can be given as below: 

hm //"" = and lim K^^"") r + irh'^'^) o] = 0. (36) 

t+r/c— const i+r/c— const 

Second, wc assume that there exists a local reference frame (with coordinates denoted as y^ = {ctbiYb)) associated 
with each body b in the A^-body system. We further assume that there exists a smooth transformation connecting 
coordinates y^ chosen in this local reference frame to coordinates x^ in the inertial frame. In other words, we assume 
that there exist both direct x^ = g^{y\^ and inverse y^ = /^(x') coordinate transformation functions between the two 
reference frames. Our objective is to construct these transformations in explicit form. 

The most general form of the post-Galilean coordinate transformations between the non-rotating reference frames 
defined by global coordinates x^ = {x^,x"'), and local coordinates y^ = (yaiUa) ^^-Y be given in the following form 

x° = y° + c-'lCa{y°, yl) + c-^Uyl yl) + 0{c~'), (37) 

^" = v: + ^?M) + ^''Qaiylya) + o{cr^). (38) 

where z^^ is the Galilean vector connecting the spatial origins of two non-rotating the frames, and we introduce the 
post-Galilean vector 2^00(2/°) connecting the origins of the two frames 

<M) = < +c-2Q^(y°,0) +0(c-4). (39) 

The functions /Ca, Ca and Q'^ are yet to be determined. 

The corresponding inverse transformations (further discussed in Sec. |V| are given by 

yl - X°+C-^]Ca{x\x')+C-^Ca{x\x')+0{c-''), (40) 

yZ = x"-zl,{x^)+c-^QZ{x\x^)+0{c-^), (41) 

where 2:^g(x°) is the Galilean position vector of the body a, expressed as a function of global time, x°. 

The local coordinates j/a are expected to remain accurate in the neighborhood of the world-line of the body being 
considered. The functions ICa,Ca and Q^ should contain the information about the specific physical properties of the 
reference frame chosen for the analysis and must depend only on the mutual dynamics between the reference frames, 
i.e., velocity, acceleration, etc. Furthermore, we would like these transformations to be smooth, which would warrant 
the existence of inverse transformations discussed in Sec. |Vl 

The coordinate transformation rules for the general coordinate transformations above arc easy to obtain and express 
in the form of the Jacobian matrix dx"^/dy". Using Eqs. ([57| - (|38p . we get: 

where Vq^ = Zq^ = cdz^^/dyl is the time-dependent velocity of the frame a relative to the barycentric inertial frame. 
We can verify that, in order for Eqs. (gOll-dll]) to be inverse to Eqs. ([371) - ([35]), the "hatted" functions (ta,Ca, Q") 
must relate to the original set of [JCa^Ca, Q") via the following expressions: 

/Ca(x'=) = -/Ca(x°,rl)+0(c-4), (44) 

QZi^') = KJc)JCa{x\rl)-Ql{x\r'^^ + 0{c-^), (45) 



with r'^a — ^'^ ~ ^a J where x^ (x") is the post-Gahlcan position vector of the body a expressed as a function of the 
global time-like coordinate x^ (as opposed to Eq. ([5^ . which is given in local time y^) defined as 

<o-<-c-'S^(^^0) + 0(c-4). (47) 

Also, w"jj = z^^ and a;'"' = (x°,a;'^). 

The transformations given by Eqs. pT ^ - pSp transform space and time coordinates from the global, inertial frame 
{x™} to space and time coordinates in the local, accelerating reference frame {y^}- These transformations are not 
singular in those regions where the post-Galilean approximation remains valid. To check this, we can verify that the 
determinant of the Jacobian matrix does not vanish [lO| : 

This guarantees the invertibility of the Jacobian matrix and hence, the existence of the inverse transformations. The 
inverse of the Jacobian matrix (|42]) - p3l) . dy^ /dx™, can be obtained directly from (|40p - (|4T|) : 

dx° ^ dx^^ feo ^ ^ '' dx<- dx'-^ ax" ^ ^ '' ^ > 

Third, we will search for the solution of the gravitational field equations, g™", in the case of A^-body problem within 
the general theory of relativity in barycentric coordinate frame in the following form: 

^ dVb dVb 
where the terms on the right-hand side have the following meaning: 

• The first term, 7™"(a;) = (1, —1, —1, —1), is the contravariant form of the Minkowski metric of flat spacetime; 



• 



The second term is a summation term representing the superposition of the unperturbed one-body solutions 
^6(0) ^^"^ each body b in the system of equations ((35)) . which were boosted to the barycentric frame by coordinate 
transformations (I37l)-(l38l): 



• Finally, the last term, /i|^", represents the gravitational interaction between the bodies in the A^-body system. It 
is expected that the interaction term is at least of the second order in the gravitational constant, /ij^" c>c G^ « c~*; 
this assumption will be used in our approximation method. 

Note that the contravariant components of the metric tensor g™" allow us to rely directly on the form of the coordinate 
transformation x™ = /™(2/a) from Eqs. (|37l) - ([38|) and the associated Jacobian matrix dx^ jdy^ given by Eqs. p2|) - 
((33]). Because of the metric nature of the general theory of relativity, Eq. ([51]) is valid for any admissible covariant 
coordinate transformations x™ = /™(j/^). As, at this stage we already know the unperturbed one-body solution /i^/g-j, 
which is given for each body h in the form of Eqs. ([35]), the only unknown part of the ansatz is the gravitational 
interaction term, /i™, which will be determined in Sec. Ill Fl 

The form of the iterative solution for A^-body metric given by Eq. (j5ip was inspired by a successful perturbation 
theory approach of high-energy physics. In that approach, a solution for a system of interacting flelds is sought in 
the form of a sum of the unperturbed solutions for each of the field in the system, plus the interaction term. The 
fact that the line element, formed using the metric tensor, represents a Lagrangian for test particles (see discussion 
in Sec. IIV Ap . allows one to develop an iterative solution in the case of A^ interacting bodies based on the principles 
of the perturbation theory approach. 

Fourth, we impose the harmonic gauge conditions given by Eq. ^, dn{^/—gg^") = 0, on the A^-body metric tensor 
girin gjygjj^ ^y gq (j5ip _ ill both the global and local frames. As we will see in the following section, this condition 
that is imposed on the metric tensor leads to the requirement that the transformation functions (|37p -(|35 )) must be 
harmonic. Below we investigate the impact of the harmonic gauge on the transformation functions ICa,Ca and Q^ 
and the inverse transformation functions /Cq, Qa and Ca- 



D. Imposing the harmonic gauge conditions 

The dynamical condition, i.e., the requirement that the spatial origin of the transformed system of coordinates 
^m _ f(^y™'^ is to move along a specific world-line, does not uniquely define y™. The existence of this coordinate 
freedom allows us to impose the harmonic gauge condition on the metric density g™" = \/—gg""^ in the local frame: 

d 



a vv-59""')=0. (52) 

The vanishing of the covariant derivative of the metric tensor (i.e., i^—gg""'^)-m = (9/i9y™)(V~55'"") + \/ ~ 99^^^1zi 



ki — 

0, where F^^ are the Christoffel symbols associated with the metric tensor g^n in the moving frame {y^}, and the 
semicolon represents covariant differentiation with respect to the coordinates, allows us to present Eq. (|52l) in the 
following equivalent form: 

g^'Tli = 0. (53) 

Furthermore, remembering the rules of transformation of the Christoffel symbols under general coordinate transfor- 
mations, we can see that Eq. (|53p is equivalent to imposing the harmonic conditions on the transformation functions 
in Eqs. dSZD-dMl): 

Uyx"" = 0, (54) 

where Dj, = {\/ — g)^^ {d / dy^){yy — gg^^'^d / dy'2) denotes the d'Alembertian with respect to {Ha} acting on x™, which 



are treated as individual scalar functions. 

Conversely, we also impose the harmonic gauge conditions on the inverse transformation: 



dx" 



35"™) =0 and 0,^ = 0, (55) 



where D^; = {^ — g)^"^ {d / dx^"){\/ — gg"^^ d / dx''^) denotes the d'Alembertian with respect to global coordinates {xm} 
acting on the scalar functions y™. 

Therefore, on the one hand the harmonic gauge imposes restrictions on the partial derivatives of the metric tensor, 
as seen in Eq. ((52|) . On the other hand, it restricts the choice of admissible coordinate transformations only to 
those that satisfy the harmonic equation (|54p . These two consequences of imposing the harmonic gauge will be used 
to establish the structure that the metric tensor, which will be given by Eqs. (|126p ~ (|128p . must satisfy under the 
harmonic coordinate transformations ([57 1) - ([55|) . The harmonic gauge also constrains the form of the transformation 
functions [JCa, Ca-, Q")- Below we explore this important dual role of the harmonic gauge in more detail. 

To study the harmonic gauge conditions we will use a complete form of the metric tensor ([T])-®, where the 

component goa has the term of c^^ order, i.e., goa = cT^UQa + c^'^Sod + C'(c~^), where ^g^ has inertial nature [16| . 
As a result, the harmonic gauge conditions ((5^ yield the following set of partial differential equations for g™": 

dylV ■ J dyl L J 

A{e-^w- + o(c-3)} + Irl^-^^'^i- + o{cr-)\ = 0, (57) 

where we use bracketed superscript indices in the form gl'^1™" to denote order terms with respect to inverse powers of 
c~^ , as discussed in the previous section. Grouping terms by order in these equations, we obtain the following set of 
equations on the metric components: 

Ajl*. = o[c-\ (58) 

where the metric density components in the expressions above can be expressed via the covariant metric as 

ffino" = -I'^'^all (61) 
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9^''"" = -I{ir7'i9fl+9N^I)}, (62) 

5'^'°" = -7"^{.S + l9'o} (7-(.f 1 + 9'^}9li) - .^J) - l''9lVA}, (63) 

^I^I"^ = -7"V^(.S-.^l'.S)+^7"^(.9?J+7^^(..S-5^t'.^^i)). (64) 
As a result, the harmonic conditions Eqs. ([58| -([60 |) take the following form: 

^ gl^} = 0{c-^), (65) 



dy, 

^c|,{5^1-7^^(^'+.M)} + 5£;{.S + |9iX^"(^-^ = Oic-'), (66) 

Ui+afW-5M-i7..(.S+7^^(5S-5^l'5S))} = 0(^-^)- (67) 



_d_ [1] _a_ 



A general solution to Eq. (|65|) may be presented in the form below: 

9l%a) = Ma(y2) + £a.(yS)ya + (^(c""), (68) 

where ^q, and £0,^ (with eav = —^va) are arbitrary functions of the time- like coordinate y°. The function /i^ may 
be interpreted as a rate of time shift in the chosen reference frame, while €ai, represents a rotation and shear of the 
coordinate axis. By an additional gauge transformation of the coordinates to a co-moving reference frame, we can 
always eliminate /i^ without loss of generality. We can also require the proper reference frame of a moving observer 
to exhibit no Newtonian rotation or shear of its coordinate axes. In addition to the fact that the leading term in 
the mixed components of the energy- momentum tensor is of the order 5*00 oc c (see Eq. (|135p below) , in terms of the 
metric tensor, these two physical conditions are equivalent to a requirement that all mixed components of the metric 
tensor of the proper reference frame at the c~^ order vanish. In other words, (/qJ, has the form: 

gW=0(c-4). (69) 

This condition leads to a significant simplification of the system of equations ([65)) -([GT ]) : 

il = 0(c-% (70) 

i4{.,s-o-.,a}+^£ = «(c-). (71) 

9 r [21 1 



dyaiL 



{.9LV^7..(..^1+7-.S)} =0(c-). (72) 



We note that, typically, researchers use only the last two of these equations. This is done to constrain the degrees of 
freedom in the Ricci tensor. For these purposes the harmonic gauge conditions reduce to the following two equations: 



iAr^- [2]_ [2]\ ^ 

2dyiy ^^'' ^°°J % 

1 d / [21 . ,„ [9i\ ^„ d 



-„{-<" 9?} -&)'^& =0(c-»), (73) 



However, as we shall see shortly, the harmonic gauge conditions offer more than is typically used. Below we shall 
explore these possibilities. For this, by formally integrating Eq. ([7^ . we get 

5iV^7"/^(5oo+7^'5S) =MoMj/") + 0(c-^), (75) 

where jiap is an arbitrary function of the time-like coordinate y^ only, representing a uniform scale expansion of the 
spatial part of the symmetric metric gmn ■ After some algebra, Eq. (j75p may be rewritten in the following form 

55 + ^o/^sS = ^^o.p{yl) - 7a^yVeA(y°) + 0{c-^). (76) 
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The asymptotic boundary conditions (|36|) require that the arbitrary function of time Hap that enters Eq. ((76)) must 
be zero. Therefore, the boundary conditions require fJ-apiVa) = Oj thus, the chosen coordinates are isotropic, we are 
led to the fohowing form of the gauge conditions ([70 |) - (f72|) : 



d 



J2] 



^"dyf^^ ■ dy. 

-,[2] 



i/Oa 
d [3] 



0{c-'), 



-,[2] 



Equivalently, for the contravariant form of the metric from Eq. (|55p . wc obtain: 

d 



o d 
^'d^' 



2100 



dx'' 



g[2]aP _|_ ^a/3^[2]00 



g[3]0. ^ 0(^-2)^ 

0(c-2). 



(77) 
(78) 

(79) 

(80) 
(81) 
(82) 



These sets of gauge conditions form the foundation of our method of constructing a proper reference frame of a moving 
observer in the gravitational A^-body problem that will be discussed next. 



E. The form of the functions of the harmonic coordinate transformations 



We now explore the alternative form of the harmonic gauge given by Eq. (|54p . Substituting the coordinate transfor- 
mations (|37| -([38 |) into Eq. (|54|) . we can see that the harmonic gauge conditions restrict the coordinate transformation 
functions ICa,Ca and Q" such that they must satisfy the following set of second order partial differential equations: 



dyf 



■7 



7 



cA 



dy^adVa 
dVadVa 



dVadVa 

Analogously, for the inverse transformation functions wc obtain: 






dx°^ 



7 



+ 7^ 



-7 



eX 



^ ^ l^g 

dx'^dx^ 
dx'^dx^ 



0{c-% 



(83) 
(84) 
(85) 

(86) 
(87) 



The general solution to these elliptic-type equations can be written in the form of a Taylor scries expansion in terms 
of irreducible Cartesian tensors, which are symmetric trace-free (STF) tensors. Furthermore, the solutions for the 
functions JCg,Cg and Q" in Eqs. (|55)) - (|85p consist of two parts: a fundamental solution of the homogeneous Laplace 
equation and a particular solution of the inhomogcneous Poisson equation (except for Eq. ()83p . which is homogeneous). 
We discard the part of the fundamental solution that has a singularity at the origin of the local coordinates, where 
y" = 0, as non-physical. 

As discussed at the beginning of this section, Eqs. ([77| -([79 |) provide valuable constraints on the form of the metric 
tensor. As a matter of fact, these equations provide two additional conditions on JCg and Q". It follows from Eqs. (|77p . 
([75)1 and the form of the metric tensor gmn given by Eqs. (|126p ~ (|128p (which are yet to be discussed in Sec. IIII A|) 
that these two functions must also satisfy the following two first order partial differential equations: 



cdy^ 



0{c- 



(89) 
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c dy'^ c dyP dy'i dy'^ \ dy'i V 

whereas the inverse functions ICa and Q" satisfy 

The two sets of partial differential equations on /Cq, Ca and Q" given by Eqs. ([83|) ~(|85 |) and (|89| -(|90 |) can be used 
to determine the general structure of JCa and Q". Similarly, Eqs. (|55)) - ([55)) and (PT|) - ([M)) determine the general 
structure of Ka and Q". These will be required in order to solve the gravitational field equations. However, we must 
first develop the solution to the A^-body problem in the quasi-inertial barycentric coordinate reference system. 

F. Solution for the gravitational interaction term the global frame 

In order to find the gravitational interaction term /i|^"(x) in Eq. (|5ip and, thus, to determine the solution for the 
A^-body problem in the barycentric frame, we will use the field equations (U) of the general theory of relativity. To do 
this, first of all we need to obtain the components of the metric tensor (|5ip in the barycentric reference system. This 
could be done by using the solution for the unperturbed one-body problem presented by Eq. (|35p . the transformation 
rules for coordinate base vectors given by Eqs.(|35])-(|33]). Using all these expressions in Eq. ([FT]) results in the following 
structure of the metric components in the global (barycentric) reference frame: 



g°°(a;) = 1 + ^ X! Y"b{yb{x)){l - -^VboAo) - -^'"boeWl{yb{x)) + 

b 

+ ^w,{y,{x))[^ + \v,,A.) + ^wl{jj,{x))] + l/iW™(x) +0(c-S), (93) 

5°"W = 4E{"'"(2/6(:^))+<o^6(y6(x-))}+0(c-5), (94) 

g-P[x) = ^-P-^'^p\Y.vo,{y,{x))+0{c~% (95) 



b 



Next, we need to specify the "source term" for the A^-body problem on the right-hand side of Eqs. (|35|) . Similarly 
to (j5ip . this term may be written in the following form: 



dx"^ dx 



^""(^) = E7rr5V^"(?/''(^))' (96) 



6=1 



dy'^ dy^ 



where 5*^' is the "source term" for a particular body h which was used to derive the one-body solution Eqs. ((T5)) - (PT]) . 
By using the source term (|96p and the transformation rules for coordinate base vectors given by Eqs. (|42|) - (|43l) . 
and definitions Eqs. (fT8 |) -(f21 ]) . we obtain the components of the source term 5™" in the barycentric reference frame: 



S''{x) = ^^^ E { (l - |-.0.<)^'>(J'^(-)) - |-.oe^6(2/.(-)) + |(|^ + ¥b.Ayb{yb{x)) + 0(c-4)}, (97) 

s'"ix) = cJ2{<^biybix)) + VboMybix)) + o{c-^)}, (98) 

b 

S-^ix) = -r^lc'Y,{a,iy,ix)) + 0{c-')}. (99) 

b 

The components of the Ricci tensor in the global frame can be directly calculated from Eqs. P^ - P^ : 

i?°°(x) = la^^cr'g^'^''^ + c-^{gi^^"'-^^{gi'^'"Y})+c-^l^^ (100) 
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i?°"(x) = c-3iA^5[3lO" + o(c-5), (101) 

R'^f^ix) = c-^^A,g^^^^f^ + Oic-^), (102) 

where A^ = j'^'^d^/dx'^dx^ is the Laplace operator with respect to {x^}. 

Substituting the components (|93| - (|95|) of the metric tensor and the components (|97| "([99 l) of the source term into 
the field equations ^ of general relativity with the Ricci tensor (|100p - (|102p , we obtain the gravitational field equations 
of the general theory of relativity in the barycentric reference frame: 

- -^\^Wb{Vb{x)yj +^/i|^l°°(x)| = 47rG^|(^l- ^i;f,o,Wbojcr6(2/6(a;))- — i;f,„,crb(2/fc(x)) + 

b b ' ' 

+ ^i^ + W,A>b{yb{x))}+Oic-% (103) 
A,^{z^,"(y,(x))+«,"„«;,(y,(a;))} = 47rG^ {<(yb(x)) + <a,(y,(.T))} + ©(c'^), (104) 

6 b 

A^Y.'^biybix)) = 47rGY,'^biybix))+Oic-^). (105) 

& b 

We observe that Eqs. (|104p and (jlOSp are satisfied identically. We will focus on Eq. ()103p with the aim to determine 
the interaction term. For these purposes, we use Eqs. (|104p and pOSp and rewrite (|103p in the following form: 



0^{j2[^biybix)} + ^{-Q4 + kvboeVlo)Mybix)) + -^wl{ybix))\--(^J2wb{ybix))'^ + 2^^!nr(a;)} = 

b ' ^b ' I, 

= 4nGj2{ab{yb{x)) + ^[^ + ^v,^Ao)<^b{yb{x))}+0{c-% (106) 



To analyze this equation, it is helpful to express the d'Alembert operator D^: via its counterpart Oy^ expressed in 
coordinates ya associated with body a. Using Eqs. (|57|) - (|55| . one easily obtains the following relation: 

dx^' '^ ^ dx^dx^ " dyf ^ ^ dy^dy^, c^ T'^" + ^ dy^Oy^ J dyH 

The terms of c^^ order on the right-hand side of this expression can be further simplified with the help of harmonic 
gauge conditions. Indeed, using Eqs. ([55)) and (|90p the expression above takes the following form: 



d^ , .A d^ &' , ,, d^ 2 fdJCa , , , X ,, d 



eA CI _ a ^x a I (OK,a 1 c \ eA C* _4 



The last equation may be given in the equivalent form: 



(108) 



°^ ^ {i + |(|| + KoXo)+0(c-^)}n..- (109) 

Therefore, taking into account that, for any given yb, the following relation holds 

ny,wb[yb[x)) = AT:Gab{yb{x)) + 0{c~^), (110) 

therefore, it follows from Eqs. (|109p and (jllOp . the following result is valid for any body h: 

n,wb{yb{x)) - A^Gab{yb{x))[l + |(|^ + !%.<,) + 0{c-^)]. (Ill) 
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With the help of Eq. (jllip we can now rewrite Eq. (|106|) as below 

°4E [|(S + Ke<)^^(2/^(-)) + i^iivbi^))] i{Y.^»{y,{x))y + ^^Tw} = o(c-')- (112) 
b ' ^b I, 

This equation can be solved for h\^^ (x). The solution to this partial differential equation is composed of two 
families: one is a particular solution that can be found by equating the entire expression under the D'Alembert 
operator to zero, h}J^(Q\{x), while the other is a family of solutions involving STF tensors that are divergent at spatial 
infinity; 

hTi^) = h\^^^^^{x)+Y^^^ShZ,,...,,{x')x^^-x'^^+0{\xT)+0{c-% (113) 

fc>0 

where ft.PJ't „ „ (a;") being arbitrary STF tensors that depend only on the time- like coordinate x". Requiring that 
the solution to this equation will be asymptotically flat and and that there is no gravitational radiation coming from 
outside the system Eq. p6p namely 

;im hTi-)-0 and ^im [{rhTlr + {rhT),o] - 0, (114) 

t+r/c— const i+r/c— const 

implies that h^^^ Uk(^*^) ^ *^' ^^^ ^^^ ^' Therefore, Eq. (|113p yields the following solution for the interaction term: 
hTi^) - 2{[j2wbiyU^))y -Y,wliy,ix))} -4j2nj,{y,{x))[^ + i«,„X„) + 0{c-'). (115) 

b b b ^b 

The first two terms on the right-hand side of the equation above represent the usual composition of the gravitational 
interaction term for A^-body system. The last term is the coupling of the local inertia field to the local gravity potential 
transformed to the global coordinates. The presence of this coupling reduces the total energy in the system, acting 
as a binding energy for the combined system of gravity and inertia. 

Substituting the interaction term Eq. (|115|) into the .g°°(a;) component (|93)l of the metric tensor we obtain the 
following expression for the metric tensor in the global frame: 

g°°(.T) = l + ^Y.nJbix) + ^(Y,Wbix)y + 0{c-'), (116) 

^ b "^6 

^ b 

2 



r^ix) 



i''''-r^-Y.'"b{x)+o{c-^), (118) 



^- b 



where the gravitational scalar, Wbix), and vector, w^(x), potentials of a body 6 transformed to the global coordinates 
{.T™} have the following form: 

Wb{x) = {l-—Vf,^y^,^)wb{yb{x))--^v^,^^wl{yb{x))+0{c-), (119) 

wt{x) = w'S{yb[x)) + v^wbivbix)) + 0(c-2). (120) 

It is convenient to express these potentials in terms of the integrals over the volumes of the bodies in the global 
coordinates. To do that we use Eqs. (jllip and (|104p to determine the following solutions for the scalar Wb{yb{x)) and 
vector w"{yi,{x)) potentials that satisfy boundary conditions (pS)) : 



j3 ,Myb{x')) , J_^£^ f ^3 , 
2c^ dxo' 



Wbiybix)) = G I d'x' "';}'''y/ +^G^ I dVa,(y,(x'))|x-x'| + 



-"/^'^'i:^(i-j%.<)-°<-''' 



(121) 



wnVbix)) = G I d3^'^£M^ + o(c-'). (122) 
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Substituting these results in Eqs. (J119p - (|120|) in the form of the integrals over the body's volume leads to the 
following expression for gravitational potentials of the body b in the global coordinates of inertial reference frame: 



».(.) . (i--|„,,.„UG/A':ifii^-l.v,,G/A'fgi^ 



C^ " / X — X' 



1 ^c-'d'' 



dVa,(,,(xO)|x-x'| + ^|dV^^;[^(|| + XXo) +0(c^^), (123) 



,3 /<(yb(a;')) , „,a^ [ ^3,^b{yb{x')) , ,^l^-2^ 



<(.t) = Gj d^x'llM^-J^+vl^Gj rf^x' "— |- +0(c-^). (124) 

We have obtained a solution for the gravitational field equations for the A^-body problem in the post-Newtonian 
approximation. In the global (barycentric) reference frame this solution is represented by the metric tensor given 
by Eqs. (|116p - (|118p together with the scalar and vector potentials given by Eqs. (|123p and (|124p . Note that the 
form of the metric (|116p - (|118p is identical to that of the isolated one-body problem ([55)) . except that these two 
structurally-identical metric tensors have different potentials. Such form-invariance of the metric tensor is due to 
harmonic gauge conditions and the covariant structure of the total metric sought in the form of Eq. (j5ip together 
with the covariantly-superimposed form of the iV-body source term Eq. (|96p . 

In the following section we transform this solution to the local frame, verify that this solution satisfies the equations 
of general relativity in that frame, and improve the form of this solution using the harmonic gauge. 

III. TRANSFORMATION TO THE LOCAL FRAME 

We begin this section by investigating the transformation properties of the covariant form of the metric tensor from 
the global frame to the local frame. 

A. Gravitational field equations in the local frame 

The metric tensor in the local frame can be obtained by applying the standard rules of tensor transformations to 
the metric tensor in the global frame given by Eqs. (J116p - (|118p : 

Ox 3x 

gmniVa) = ^-^-Q-^ gkl{x{ya))- (125) 

Using the coordinate transformations given by Eqs. ([57 )) - ([55|) together with Eqs. P^ -iP^ ]) . we determine the form 
of the metric tensor in the local reference frame associated with a body a: 

2 XT- r/ 2 . \ , , .. 4 



XI { (l ^ -^^aoeV'aa)wb{x{ya)) + -^V^^^wl{x{ya))^ 
b 



+ |(E^^(^(y")) -^- \v^,,.vl)\o{c-% (126) 



The form of the metric in the local reference frame given by Eqs. (|126p - (|128p still has a significant number of degrees 
of freedom reflecting the presence of the yet arbitrary transformation functions. Some of these degrees of freedom 
can be eliminated by imposing the harmonic gauge conditions. Eqs. (j77p and (|79p can provide valuable constraints 
on the form of the metric tensor in the local frame. Indeed, using the actual structure of the metric expressed by 
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Eqs. (|126p - (|128p . the gauge conditions (|77l) - ([79|) and Eqs. ((89)) - (|90p . we can now use these two sets of constraints on 
the transformation functions and to improve the form of the metric tensor in the local frame. As a result, the metric 
tensor in the harmonic coordinates of the local frame has the following form: 






+ ^,{T.^'>«ya))-^-h-ao.<J + oic'% (129) 



b 

^ b 

1A 



gafiiVa) = lafi+lafi-^y^Wb{x(ya))- ^-^ -\v^^y^^^+0{c^'^). (131) 



b 

The contravariant components of the source term in the local reference frame, SmniVa), are related to that in the 
global frame, Smni^), via the usual tensor transformation: 

dx dx 

SmniVa) = Q^ Q^ Skl{xiya)). (132) 

On the other hand, given the definition for the scalar and vector densities (which are given by Eqs. ((T^ - (PT|) ). for 
which we need covariant components of the source term S"™", we can determine Smniya) in the local frame via the 
covariant components 5*™" in the following form 

Smn{ya) = 9-nik{ya)gnl{ya)S^\ya) = g,nk{ya)gnl{ya)-g-J -g^ SP''{x{ya)), (133) 

where dy^/dx^ are given by Eqs. (|49|) ~([50 | . which lead to the following expression for the components of the source 

term Smniya)'- 



SooiVa) = {l + ^gl^{ya) + 0{c-'))x 

^ 2^'^ {(^ " ^«aoe<„)fTfc(x(ya)) + -^Vaoealixiya)) + ^(^ + 5^a„e<,) ^^(^^(ya)) + 0(c-^) } , (134) 
b 

Soo.{ya) = j^^cJ2{<^b{^iya))-via,{x{ya)) + 0{c'')}, (135) 

b 

So^piVa) = ~^ap^c'[Y,<^b{x{ya))+0{cr')], (136) 

b 

where ICa comes as a part of the rules Eqs. (|49| -(|50 l) that are developed for the inverse transformations j/™ = f™{x^). 
Similarly, to compute the Ricci tensor in the global frame, we use the following expression: 

Rmn{ya) = gmk{ya)9nl{ya)R^\ya), (137) 

with R^^{ya) being expressed in terms of the covariant metric gmn, R^''[gmn{ya)], namely: 
RooiVa) = [l + ^gl^^ix)+Oic-')) 

X ( - In,^ (c-.gi^J + c-{^i - i(,t^i)^}) + cr^lg^/o^ {^1 + 7..5S) + Oic-^)) , (138) 

RoaiVa) = -c-'^AyJ^i+Oic-'), (139) 



RMVa) = -c-^lAyJ^l + 0{c-^), 
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(140) 



where Aj^^ = l'^ d /dyl^dy^ is the Laplace operator with respect to the coordinates {y^}. 

The metric tensor given by Eqs. (|129p - ()13ip and the condition given by Eq. (|79p allow us to present the gravitational 
field equations @ of the general theory of relativity in local coordinates in the following form: 



^v\H{[^-^''a,c<,)M<ya)) + -^^ao,<i<ya))} 



V„„^V„ 



2 



2 fdCa , i/a/Ca\2 , OQl fdICa , i 






dyl) ^"^dy^, KdyO 

4 _ , .. 2 ,dlCn 



A„ 



47rG5^{(l--.;,„Xo)^fc(^(ya)) + ^^ao.^^Wya)) + ^(^ + i^'.^.^Ja.Cxly^)) + ©(c-^)}, (141) 

b 



cdy- ' ^-"92/2 ' ""'dy- "°" 9y0 



A.„[E^''(^(y«))- a^ - \va,.<^+0{c-^)\^S^G^^pY.^,{x{ya))+0{c- 



dyl 



(142) 
(143) 



Using the fact that A^^ == A^; + 0{c~^) and the gauge equations (l83])-(l85]) we can show that Eqs. p42p and p43p 
are identically satisfied. Furthermore, Eq. p4ip reduces to 

□y. E^''(-^(y-)) = 4^GE {^^(^(ya)) + |(|^ + 5«ao.O^(^(ya)) + 0{c-*)}. (144) 

b b 

We can verify the validity of this equation in a manner similar to the derivation shown by Eqs. p07p and (jlOSp : 
We can express the d'Alembertian in the local coordinates Dj,^ via global coordinates {a;™} as below: 



d^ 



+ 7 



£A 



52 



52 



„eA 



92 



dyOj dyldy^ dx^' 



^{"^ 



M -y^ 



92Q^ I d 



}; 






,9Q^ .,9Qn 92 



+ r 



+ 0{c-'), (145) 



where dy^/dx'^ needed to derive Eq. (|145p are given by Eqs. P^ - ((5(I)) . Terms of order c^^ t]^at arc present on the 
right-hand side of this equation can be simplified using the harmonic gauge conditions ([55]) and (IMI) , so that Eq. (|145p 
takes the form: 



92 



+ 7 



92 



92 



ayf a^eg^^A ^^0 



+y 



92 



2 /a/Ca , 1 



^'aoe<o) 






The last equation may be given in the equivalent form: 



(146) 



(147) 



As a side result, comparing this expression with Eq. p09p . we see that the following approximate relation holds: 

(S + ^-^oe<) --{^ + '.-..<) + ^(^-^)- (148) 



Using Eq. (|3T|) we see that the following relation satisfied for any {ya}'- 

a^Wbixiya)) = 47rGafc(.T(y,))+0(c-4), 



(149) 



18 
and relying on (|147|) . wc observe that the foUowing relation also holds: 

Uy^w,{x{y,)) = A^Ga,{x{ya))[l + ^[^ + \v,^,vl^)+0{c-^)], (150) 

which proves Eq. pli)) . As a resuh, with the help of Eqs. (IH31)-(I55]), and also ([T50| . we verified that Eqs. ([TiT|) - pi5)) 
are satisfied identically. 

In the next subsection, we use the harmonic gauge conditions to constrain this remaining freedom which would lead 
to a particular structure of both the metric tensor and the coordinate transformations. 

B. The form of the metric tensor in the local frame 

Eqs. (|129p - (jl3ip represent the metric tensor r]mn given by Eqs. ([771) - ([75|) in the coordinates of a local reference frame 
that satisfies the harmonic gauge. This set of gauge conditions forms the foundation of our method of constructing a 
proper reference frame of a gravitationally accelerated observer. As a result, the metric representing spacetime in the 
proper frame may be presented in the following elegant isotropic form that depends only on two harmonic potentials^: 

gmiVa) = l-\w{ya) + \w-'{ya) + 0{c-^), (151) 



4 



50a(2/a) = -7aA3U'^(2/a) + 0(c-5), (152) 

2 






gapiya) = laP+laP^w{ya) + 0{c ), (153) 

where the total (gravitation plus inertia) local scalar w{y^) and vector w°'{y^) potentials have the following form: 

W{ya) = 2^Mya)-^-kv,^Xo-^{-Q^+^^aoe^ + l[-Q^) -(^ + 5«ao.<oj | + ^(c ^),(154) 

-"(..) - j:<iy^)-l{rf^^^^'^^r^^.J^-<'^}^oic-^), (155) 

with gravitational potentials Wb{ya) and w'^{ya) transform between the global and local frames as below 

Wb{ya) = [l- -^v^^y^^^wb{x{ya)) + -^v^^^wl{x{ya)) + 0{c^^), (156) 

<{ya) = w<^{x{ya))^vlwb{x{ya)) + 0{c-^). (157) 

It follows from Eq. (j78p that potentials Wb{ya) and w'^{ya) satisfy the continuity equation: 

c^^^'-^^Oic-^). (158) 

dyl dy^ 

With the help of Eqs. ([55|) - ([55|) and ((M)) ~ (pn)) . we can verify that the new potentials w and w" satisfy the following 
post-Newtonian harmonic equations: 

^VaWiya) = 47rG^ |(^1 - -^VaaeVl^jcTbixiya)) + -^V ^^^al{x{y a)) + -^{-Q^ + \'^'aaeV'a„)(^b{x{ya)) + Oic-"^)^, 

(159) 
^y^w-^i^ya) = 47rG^ {a,"(a;(ya)) - vl^ab{x{ya)) + 0{c-^)]. (160) 



Alternatively, one can write the metric tensor, given by Eqs. I ll51l |- ||153t . in the following equivalent harmonic form: 
1 + c~'^w{ya) c-^ \ / 
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One may verify that according to Eqs. (|150p and (J142I) . the scalar and vector potentials can be presented as functions 
of the local coordinates only: 



.2a2 






2G f,,,a,ix{y^fdlCa^, 



c2 






«;aa;(ya)) = G / d^y^^^^i^ + 0(c-^). (162) 

These expressions allow us to present the scalar and vector potentials, given by Eqs. p56p - (|157p . via barycentric 
densities given as functions of the local coordinates: 

,<^b{x{y'a)) 4 /"a ,<Jl{x{y',)) 






,3„/ <(a;(ya)) „,a^ [ .fi„.l^b{x{y'a)) , ^^^„_2^ 



«;,"(y.) = cy ^'^i |;;_;;i -<ogy ^ ^^f"^"''^^^^- ^'^^^ 

One can also present the transformation rules for the relativistic gravitational potentials from the proper coordinates 
{yb] associated with body h to the coordinates {ya] associated with body a. With the help of Eqs. (|119p and (|120l) . 
we can express the potentials Wb{x) and w^{x) as functions of {x{ya)} and substitute the results into Eqs. p56p and 
P57p . as shown below: 



Wb{ya) = [I- -^Vba„e^tao)wb{yb{ya)) + -^Vba„^wl{yb{ya)) + 0{c ^), (165) 



2 \ 4 

^Vba,Aao)Myb{ya)) + ^ 

wt{ya) = wUyb{ya))~vt,^^Wb{yb{ya)) + 0{c-^), (166) 



with v^^ = v"^ — w" being the relative velocity between the two bodies. Assuming the existence of the coordinate 
transformations between the two bodies in the form 

y° = y", + cr'lCba{y"a,y:)+0{c~^), (167) 

Vb = y: + z,",^(y°)+0(c-2), (168) 

with function JCba given Q as ICbaiya^yl) = ^a{y^,ya) — ^b{y^,ya + zl^^) + 0{c^'^) , we note that the scalar and vector 
potentials, Wbiybiya)) and w^{yb{ya)) correspondingly, satisfy the following equations: 

Oy^Wbiybiya)) = ^7rGab{yb{ya)){l ~ |(^ + k^baoAa.) + Oirr*)}, (169) 

Ay^w^{yb{ya)) = A7rG<j^{yb{ya)) + 0{c-'). (170) 

In Eq. (|169p . similarly to Eq. (|148p we used the following identity: 

Eqs. (|169p and (|170p have the following solutions 

WbiybiVa)) = ^/^'^^ "[yj-^l + ^G^ J d^y',ab{yb{y'a))\ya y'a\ 

w^ybiya)) = G fdy ff^y'}y'^U o{c-^). (173) 
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Therefore, finally, we obtain the following expressions for the scalar and vector gravitational potentials, expressed 
as integrals over the compact volumes of the bodies given in the local coordinates: 



^ J W a 



y,| ■^^-o^^i^^nya-y; 



2Q2 



d^ya'^b{yb{ya))\ya-ya\--J- / '^ J/aT— — -7t(^0-+ 2«6aoe"fcao) +C'(C ■'),(174) 



<iya) = cfd^y'f-fMim-v^^^Gfd^y'/J^^ (175) 

We use the phrase harmonic metric tensor to describe the metric tensor p5ip - ()153p . defined in terms of the 
harmonic potentials w and w"', which in turn are given by Eqs. (|154p - (|155p and satisfy the harmonic equations 

C. Applying the harmonic gauge conditions to reconstruct the transformation functions 

The second role of the harmonic conditions in our derivation is to put further constraints on coordinate transfor- 
mations, taking us closer to determine in full the form of the functions /Cq, Q" and Ca- 

1. Determining ttie structure of Ka 

The general solution to Eq. (|83p with regular behavior on the world- line (i.e., omitting terms divergent when 
ly^jl — > or solutions not differentiablc at |y^| = 0) can be given in the following form: 

lCa{ya)^Ka,+Ka,^y':,+5Ka{ya)+0{c-^), where ^^a (^a) = ^ TJ^aoMi ...M. (2/a)yf -yf + C^(c"^) , (176) 



k=2 



with Kaof_i-i^...fi^.{ya) being STF tensors 2J|, which depend only on the time-like coordinate yJJ. Substituting this form 
of the function ICa into Eq. (|89p . we the find solutions for Kaofi and Kaiii...fj.k- 

«aoM = -C"aoM+^(c"')' «a„pi...M, =0(C-^), k > 2. (177) 

As a result, the function ICa that satisfies the harmonic gauge conditions is determined to be 

This expression completely fixes the spatial dependence of the function ICa, but still has an arbitrary dependence on 
the time- like coordinate via the function Kao(ya)- 

2. Determining the structure of Q^ 

The general solution for the function Q" that satisfies Eq. ((85)) may be presented as a sum of a solution of the 
inhomogeneous Poisson equation and a solution of the homogeneous Laplace equation. In particular, solutions with 
regular behavior in the vicinity of the world-line may be given in the following form: 

Q^iya) = < + <X + baof^r^yiiya + ^^iya) + 0(c^'), (179) 

where q" can be determined directly from Eq. (|55|) and the function 6^" satisfies Laplace's equation 

We can see that Eq. (|85p can be used to determine 9a„„i,, but would leave the other terms in the equation unspecified. 
To determine these terms, we use Eq. (^0)) together with Eq. ([M]) . and get: 

BO^ BO^ / BK \ 

^ao.-ao;3+7.A^+7^A^ + 27.,(^ + K,Xo) = 0{c-% (181) 
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Using the intermediate solution (|178p for the function JCa in Eq. (|181|) . we obtain the foUowing equation for Q": 

rlO^ BO^ /Bk \ 

z;a„a«a„,+7aA^+7/iA^ + 27.,(^ + i«,„,<,-a,„,y:) = 0{c~\ (182) 

A trial solution to Eq. p82p may be given in the following general from: 

+ y,,< + 5C(2/a) + 0(c-2), (183) 

where q^ and the antisymmetric matrix w"*^ — — w^" are functions of the time-like coordinate y^; Ci,...,C5 are 
constant coefficients; and 6^1^ {y a), given by Eq. (|180p . is at least of the third order in spatial coordinates y^, namely 
^iaiVa) oc 0{\y^\^). Direct substitution of Eq. p83p into Eq. (|182p results in the following unique solution for the 
constant coefficients: 

ci=--, C2 = 0, C3 = --, 04:^1, C5==-l. (184) 

As a result, the function Q" has the following structure: 

Q^aiVa) - <Za„ {Ho^o + < + ^"(^ + l^aoxO)yae + «ao.(2/a"y: " h^'VaXVa) + SCiVa) + 0{c-'), (185) 

where q^ and uj^^ are yet to be determined. 

By substituting Eq. (|185p into Eq. (|182p . we see that the function S£^"{ya) in Eq. (|185p must satisfy the equation: 

We keep in mind that the function S^^{ya) must also satisfy Eq. (|180p . The solution to the partial differential equation 
(|180p with regular behavior on the world-line (i.e., when |y^J — > 0) can be given in powers of y^ as 

5C{ya) = Y.h^^ao,^...,M)y'a'■■■y'a'+0{\yii\'')+0{c-'), (187) 



k>3 



where <5CaoMi---Mfc(2/a) ^^'^ ^TF tensors that depend only on time-like coordinate. Using the solution (|187p in Eq. (|186p . 
we can see that S^q is also antisymmetric with respect to the index a and any of the spatial indices fii...fik- 

Combination of these two conditions suggests that S(,q Uk ^ ^ '^'^^ ^^^ k > 3, thus 



6aiya)^0{c-'). (188) 

Therefore, application of the harmonic gauge conditions leads to the following structure of the function Q": 

where q" , uj^^ and Kao &re yet to be determined. 



QaiVa) = qZ - {Ho< + < + 7"1§^ + -2VaoxO)yae + ^aoe {v^Vl - b'^'yaXVa) + 0{c-\ (189) 



3. Determining the structure of La 

We now turn our attention to the second gauge condition on the temporal coordinate transformation, Eq. 
Using the intermediate solution (jl78p of the function /Ca, we obtain the following equation for La'. 

^''ife = -^^+o(=-^)-("«.<.+<.»:)-=^4(^ + K....:.)+°(=-^ (««) 

The general solution of Eq. (|190p for La may be presented as a sum of a solution 8 La for the inhomogeneous Poisson 
equation and a solution (5£aO of the homogeneous Laplace equation. A trial solution of the inhomogeneous equation 
to this equation, 8 La-, is sought in the following form: 

bLa = ck^{Va,,al,){VaX) + ^^2{<eyl){y a.yl) ~ ^3c'^(^ + i^ao.^o) (y-^a) + ^(c"')' (191) 
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where ki,k2, k^ are some constants. Direct substitution of Eq. (|19ip into Eq. (I190p results in the following values for 
these constants: 

fci = i k2 = ^, fc3 = i. (192) 

6 10 6 

As a result, the solution for 5Ca that satisfies the harmonic gauge conditions has the following form: 

The solution for the homogeneous equation with regular behavior on the world-line (i.e., when |y^| — > 0) may be 
presented as follows: 

Caoiva) - 4„(2/°) + 4oa(2;°) yl + \L.x^.{yl) vlv'a + SUva) + 0{c-% (194) 

where f aoAp is an STF tensor of second rank and 5(^a is a function formed from similar STF tensors of higher order: 

fe>3 

Finally, the general solution of the Eq. (|190p may be presented as a sum of the special solution 6Ca of the inhomo- 
geneous equation and the solution Cao of the homogeneous equation A£a = 0. Therefore, the general solution for the 
harmonic gauge equations for the function Ca{ya) = C,aa + 5Ca has the following form: 



Ca{ya) = U+Uxya+¥a.X^.yay'i+5ta{ya) + 

dy°a ^ dyl 



+ H^ao^OiVa^yD + M<.yl){ya.y:) - ¥'^,{^ + ba^.^^Yvauy:) + 0{C~^). (196) 



D. Harmonic functions of the coordinate transformation to a local frame 

As a result of imposing the harmonic gauge conditions, we were able to determine the structure of the transformation 
functions /Ca, Q" ^^id Ca that satisfy the harmonic gauge: 

ICaiya) = «ao-c(«,„X)+'^(c"')' (19^) 

QZiya) = Co - (l«o + < + T"' (^ + baoX<))ya. + «ao. (s/a"?/: " h"'yaA2/a) + (^(c"'), (198) 
/:a(ya) = 4o + 4oA y^ + \ia,X^. Vay^ + Sia{ya) + 

+ ^KkoeO - -^^(^ + ^-ao.^o)) (y..2/:) + T-0«^yl)iya.y:) + 0{c-% (199) 

Note that the harmonic gauge conditions allow us to reconstruct the structure of these functions with respect to 
spatial coordinate y^. However, the time-dependent functions Koo , 9"^ , w^^J^ , (-ao,(-ao\,(-ao\ii, and (J^aoMi-Mfc cannot 
be determined from the gauge conditions alone. We need to apply another set of conditions that would dynamically 
define the proper reference frame of a moving observer, thereby fixing these time-dependent functions. This procedure 
will be discussed in the following section. 

IV. DYNAMICAL CONDITIONS FOR A PROPER REFERENCE FRAME 

As we saw in the preceding section, introduction of the harmonic gauge allows one to determine the harmonic 
structure of the transformation functions /Co, 'Co and Q", but it was insufficient to completely determine these 
functions. To proceed further, we investigate the dynamics of a co-moving observer in the accelerated reference frame. 

An observer that remains at rest with respect to the accelerating frame does so because of the balance between 
an external force and the fictitious frame-reaction force. If the external force is universal (in the sense of the weak 
equivalence principle) and in the case of complete balance of this fictitious frame-reaction force, there will be no net 
force acting on the observer allowing it to be at rest in what we shall call its proper reference frame. The motion of 
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the observer in this frame, then, wiU resemble a free-fah that follows a geodesic with respect to the metric gmn- It 
would, thus, be natural to require that in order for the observer to be at rest in its proper reference frame associated 
with the local metric, the observer's ordinary relativistic linear three-momentum should be conserved. We can explore 
these conditions by writing down the Lagrangian of a test particle that represents the observer and find that it is 
possible to eliminate all the remaining unknown components oi K-a, Ca and Q^. 

A. 'Good' proper reference frame and local metric conditions 

The metric tensor gmn, given by Eqs. (J151|) - p53p allows one to study the dynamics of the reference frame that 
moves in response to the presence of the external force. The Lagran gian L of a test particle that moves in this system 
can be obtained directly from the metric gmn and written as below |16l . |20| : 



= -mc^jl + c-^Uv^v' -w\ + c-'^Uw^ - livev'f - 4v,w' + ^v^v'w + 0{c-^)y (200) 

The Lagrangian given by Eq. pOO|) leads to the following equation of motion: 
c-^{w" (l + c-\3w - ^v,v') - ic-^w"' + 0{c-^)) } = -d'^wU - c~2(|t)ew' + w)} - ^ved^w' + ©(c"'^) .(201) 

It can be shown that Eq. (|20ip is equivalent to the geodesic equation for the local reference frame that is freely 
falling in the external gravitational field. However, Eq. (|20ip has the advantage as it allows to separate relativistic 
quantities and to study the motion of the system in a more straightforward way. 

The condition that the test particle is to remain at rest with respect to the accelerating frame, then, amounts to 
demanding conservation of its three-dimensional momentum p" = dLtat/dv": 

p" = w" (1 + c'^(3w - ^v,v')) - Ac.-^w"' + 0{c-^). (202) 

In other words, the total time derivative of p" must vanish. We assume that the observer is located at the spatial 
origin, y" = 0. Then this condition amounts to requiring that the first spatial derivatives dfjWb and d^w^ are equal 
to that produced by the gravity of body a only, gmniVa)- Thus, we require the following two conditions to hold: 

hm -^zz;(yj = —^Wa +0{c-^), (203) 

2/"^o dyP dyZ y-=o 

hm^— -«;"(y,) = —-w^ +0(0-'). (204) 



^'"" P,„fi 



y^^o dy'l dy'l 



yo=o 



We also require that on the world-line of the body its momentum (|202p is to depend only on the potentials produced 
by the body itself and is independent on external gravitational fields, namely: 

P"\ya = w" (1 + c-'(3u;a - \v,v')) - Ac~^wZ + ©(c"*). (205) 

Thus, we require that another pair of conditions hold: 

lim «;(y,) = wA +0{c~^), (206) 

ya— >0 ^ya—"^ 

lim«;"(y,) = <| „ -t- ©(c-^). (207) 

Physically, we choose the coordinate system {y™} in the proper reference frame of the body a such that along the 
world-line of the body the external gravitational potentials in Wb and w", b ^ a and the first spatial derivatives dpwb 
and dpw^ are absent, so that the spacetime of the A^-body system in the local coordinates along the world-line gmn{ya) 
reduces to the spacetime produced by the gravity of body a only, gmniVa)- To put it another way, on the world-line 
of the body a, the presence of the external gravity will be completely compensated by the frame-reaction potentials. 
Consequently, in accordance with the Equivalence Principle, the external gravitational field in the body-centric metric 
will be manifested only via relativistic tidal effects in the vicinity of the world- line of body a. These requirements 
constitute a generalized set of Fermi-conditions in the case of the A^-body system. 
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Therefore, summarizing Eqs. (|203p ~ (|204p and p06p - (|207p . we require that the foUowing relations involving scalar 
and vector potentials w and w" hold along the world-line: 

lim wiva) = Wa\y^=o + Oic-''), lim —-w{ya)^—^Wa +0(c-4), (208) 

lim w'^iy,) = <|,„=o + 0(c-'), lim -^w^{y,) = ^w^ +0{c-'-). (209) 

lyaKo laaKo dy'a dya V'^=^ 

Applying these conditions yields additional equations that allow us to determine uniquely the explicit form of 
the coordinate transformation functions ICa,^a and Q", allowing us to fix the time in the local coordinates. As 
it was demonstrated in [ifll, the conditions (|208p - (|209p can be used to determine frame-reaction potentials. These 
potentials are the combination of the terms on the right-hand sides of Eqs. (|154p - (|155p that depend only on the 
coordinate transformation functions ICa,Ca and Q" and their derivatives. In addition, this procedure also enables 
us to determine the equations of motion of an observer with respect to the inertial reference system (logical, but 
somewhat unexpected result), written in his own coordinates {y^} = iVaiya)- These equations of motion will provide 
the dynamical conditions necessary to fix the remaining freedom in the coordinate transformation functions ICa , Ca 
and Q'^. 

B. Application of the dynamical conditions to fix the proper reference frame 

Imposing the generalized Fermi conditions (j208p - (|209p on the potentials w and w", which are given by Eqs. p54p - 
P55p . and denoting u to be the value of any function u on the world-line of body a via 

lim u{ya) = u|y„=o = "", (210) 

results in the following set of partial differential equations set on the world-line of the body a: 

E-^.-^-K„<.-^{§| + c„„,^H-i(^)^(^ + !.„,<.)'} ^ Oi.-'), (211, 

sr^dwb I f d^Cg d^Q% (dna,, ^ \\ .^. -4n /^ioa 

l^Trj + ''a,p--a[Trj^^+c^aoe-^— + a^J-—^ + Vg^^,Vg\\ = 0{c ), (212) 

b^a^V'^ c ^dyildyO dy^dy^ ^ oVa '> 



Equations (|21ip - (|214p may be used to determine uniquely the transformation functions ICa,Ca and Q". Indeed, 
from the first two equations above, (|21ip and (|212p . wc immediately have: 

Et^^-^TT-KoXo =Oic-% (215) 

where we, for convenience, have split the scalar potential Wb into Newtonian {Ub) and post-Newtonian (Swb) parts: 

Wb^Ub + \6wb + 0{c-^). (217) 

This way, that gravitational scalar and vector potentials in the proper reference frame of the body a defined by 
Eqs. (|218p - (|219p have the following form as functions of local coordinates: 

w,{ya) = Ub{x{ya)) + ^ [Swbixiya)) - 2{v,^,vl^)Ub{x{ya)) + ^v,^,wl{x{ya))) + 0{c-'), (218) 
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w'^iVa) = <(x(y,))-<„(7fc(x(y,)) + 0(c-2). (219) 

Equations (|178p and (j215p allow us to determine k^^ and present the solution for the function Ka as: 



^aiVa) = [\T.^b- h^aoAo)dya " civ^o^a) + 0{c-% (220) 

From Eq. (|216p we also have the wcU-known Newtonian equation of motion of the body a: 

riTj 

<, = -l^'Y.^ + 0{c-% (221) 

Substituting the expression (|215p into Eq. p98p . we can determine the function Q": 

QZiVa) - <, - [h« + 7"^ E ^^ + <)2^- + «ao.(2/ay: " h'^'VaXVl) + Olc"'). (222) 

b^a 

Finally, the general solution for the function Ca, given by Eq. (|196p . now takes the following form: 

/:a(ya)=4o+4„Ay.' + i4oA^yay^-^c(Ec|^-^ao.<o)(ya.y:) + lIjC(<,y:)(y,X) + '^4(ya)+0 (223) 

b=/=a ^^ 

The next task is to find the remaining undetermined time-dependent functions present in Q" and £q, as given by 
Eqs. (|222p and (|223p . To do this, we rewrite the remaining parts of Eqs. (|21ip - (|214p as a system of partial differential 
equations, again set on the world-line of the body a: 



b^a 



QyO aoe QyO ' n 9y0 y V Sj/g 



S5Z*'='-{s7^ + -«.SS:7: + ''-(lf ^""'"")} - <^(''"'' (^^^' 






,,/^'H^ -a^ + ^a^+7 -aoA^-a.^1 = (^^(c ), (226) 

^ dy^ ^cdy^dy'a oy'^oyS. oy^dy'a > 

These equations may be used to determine the remaining unknown time-dependent functions iam ^oqA, ^aoXin ^nd 
also q'^^ and i^qq , which are still present in the transformation functions. Thus, substituting the solutions for ICa and 
Q", given respectively by Eqs. (j220p and ([222]), into Eq. (|224p leads to the following solution for lao'- 






b^a b^a b^a b^a 



Next, Eq. f|225p results in the following equation for i'^ : 
¥ao-Y.{d^^^^+^^a^^§^)~<T.^''+^(^^^^^^^^^ (229) 

^T^a ^T^a b^a 

From Eq. (|226p we can determine i^^ : 

-Xo = -C+4E*?-2<,Ef^^-«aoe<+0(c-'). (230) 

b^a b^a 

Eq. ((227)) leads to the following solution for ^^f : 

b^a^y^f^ b^a ^^^ 
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The quantity £"^ is an STF tensor. The expression on the right-hand side must, therefore, be also symmetric. This 
can be achieved by choosing the anti-symmetric tensor uj"l^ appropriately. This can be done uniquely, resulting in 

<f - -V« - «.) + 21: (g - g) + Oi,-% (232, 



which leads to the following solution for i"^: 

¥t = 2 E (1^ + 1^) + 17"^ E # + (|7"^-.„.<„ + <,«fo + «) + 0{c-\ (233) 



Also, using Eq. (|232p . one can present a solution for ^"^j, Eq. (|229p . in the form 



Finally, Eqs. (|229p and (|230p allow us to determine the equation for g"^. Indeed, differentiating Eq. (|230p with 
respect to time and subtracting the result from Eq. (|229p . we obtain the following equation for g" : 



{Ho< 



- w„ 



«ao.-KoKe<o)+0(c"')- (235) 






Therefore, we can use Eq. (|235p to completely determine the function Qq. 

The true position of a body a includes terms to all orders, not just the first-order (Galilean) term. This led us to 
introduce the post-Galilean vector, x" (y^), defined by Eq. ([55)1 . Combining this definition with Eqs. (|216p and (|235p . 
we can now relate the magnitude of the frame- reaction force (acting on a unit mass) , written in the local reference 
frame of the body a, to the acceleration produced by the external gravitational field x"^ = z"^^ + c^'^Qag + 0{c^^): 

+ ^(^i:# - Ko.E ^ - K. i:#) + «(-')^ (^^e) 

"^ b^a ^y^ b^a^y^'^ b^a "^^^ 

The equation of motion (|236p establishes the correspondence between the external gravitational field and the 
fictitious frame-reaction acceleration i" {]f) that is needed to keep the body at rest in its proper reference frame. 
This frame-reaction force balances the effect of the external gravity acting on the observer co-moving with the body. 

C. Summary of results for the direct transformations 

In the preceding sections, we sought to write the general post-Galilean transformation from a global, incrtial frame 
of reference to a local, accelerating frame in the form 



xO = 2/° + c-2/C,(y,^) + c-4£a(2/^) + 0(c-«), (237) 

^" = 2/a" + <o(ya) + C-'Q:(ya)+0(c-')- (238) 

To determine the unknown functions /Ca, Q", and £a, we used the following approach: i) we imposed the harmonic 
gauge conditions on the metric tensors in the global and local reference frames; ii) we ensured that the non-inertial local 
frame is non-rotating; iii) we introduced a co-moving accelerating observer at rest with respect to its proper accelerating 
frame of a body a; iv) wc required that a co-moving observers's ordinary three-dimensional linear momentum be 
conserved on the world-line occupied by the accelerating frame. 

Together, these conditions were sufficient to determine /Cq, Q", and Ca unambiguously. /Cq and Q" are given by: 

K.a{ya) = y^(E^^-H„Xo)rf2/;°-c(«,,,y:)+o(c-4), (239) 

^OQ b^a 
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Q^aiVa) = qZ - {h« + 7"^ T.^b+ <)yae + ^^^.{y^yl - h^'yaXVa) + 0{c-^), (240) 

with the anti-symmetric relativistic precession matrix u^^^ given by Eq. (j232p and the post-Newtonian component of 
the spatial coordinate in the local frame, q^^^, given by Eq. (|235p . 

To present the last transformation function, Ca{ya), we use the explicit dependence of the Newtonian and post- 
Newtonian components of scalar and vector potentials on the body-centric potentials, as given by Eqs. p56p and 
(|157p . As a result, CaiVa) is given by: 

Mo 
where the functions £aoj ^"o ^^^'^ ^af ^.re defined using the following expressions: 

b^a b^a b^a b^a 

¥ao = -< + 4 E < - 2<, E ^^ - ^^0. < + ^('^"')' (243) 

b^a b^a 

¥f. = 2 E (1^ + 1^) + h-' E # + {K'^a..<^ + <A. + .la^„) + Oic-). (244) 

Substituting these solutions for the functions ICa , Q" and Ca into the expressions for the potentials w and w°' given 
by Eqs. (|154p - (|155l) . we find the following form for these potentials: 

w{ya) = ^wb{ya)-^{wb + yl-^ 

b b^a *^° 



[kyaya [l<iXaaof,a^„ + aao."aoA + "^^oe^aoX + '^^aoe^oX + 7eA E ^^ 



+ 2 E (^ + fr)] + U<ya)iya,y!i) + OoSla} + 0{c-% (245) 

W"{ya) - E<(y«)-E(<+2^a^)-lIl{3ya2/:-7"^yaX}«aoe-i^7T^^4+C^(c-')- (246) 

b b^a "y^^ "y^'"' 

Substitution of these solutions for the potentials into the expressions for the metric tensor given by Eqs. (|15ip - (|153p 
leads to the following form of the local metric tensor: 



b b^a ^° b b^a ^'^ 



+ \ylya \le\aa„^a'^„ + aaoeflaoA + 2aaoe"aoA + SUaoeOaoA + 7£A ^Ub + 



b^a 

dyl dy^ 



+ 2E(^ + ^)] +lI)«e2/a)(y.X)+5o%}+0(c-«), (247) 

b^a 
gociya) = ~l.X^{Y.^^(y-)'T.i'^i+ya^)"U^yaya~^''yaX}<e} + \j^S£a+0{c-')l24^^ 

^ b b^a y^ ' c uy^ 

gMya) = lc.fi+l<.0^{Y,Wb{ya)-Y.{'^b + yl^)}+O{c~^). (249) 

"^ b b^a y°- 

All terms in this metric are determined except for the function S£a, which remains unknown. We note that the 
potentials Eqs. (|245p - (|246p depend on partial derivatives of Sia- The same partial derivatives appear in the temporal 
and mixed components of the metric (|247p - (|248p . The presence of these terms in the metric amounts to adding a 
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full time derivative to the Lagrangian that describes the system of the moving observer. Indeed, separating the terms 
that depend on S£a from the Lagrangian eonstructed with Eqs. (|247p - (|249p . we get: 



2 r dSi, v^ dSi, . 2 dS£, 



(250) 



As a result, the remainder of the gauge transformation depending on S£a will not change the dynamics in the system 
and, thus, it can be omitted and, so that the potentials w and w" take the form: 



KVa) == '^Wb{ya)-'^[wb+yl 



b^a 



dwb 



^{Hya[ 



leXa-aoA + "ao.«aoA + "^^o^^aoX + "^^aaeKoX + 7eA ^ Ub 



b^a 



■n 

b^a 



dwhx dw. 



'bX , UWbe 



dyl dy^ 



+ U<ey'a){yaX)}+^^<'~')^ 



dw? 



\ya) = E<(y")-E(<+2^a^ 



T-A^v^vl - r^yaXX^ + ^(c"')- 



b b^a 

We can now present the local metric in the following final form: 



500 



iv^) - i-|{E-^(^^)-E(^^+^:|^)}4{[E^^(^^)-E(^^+^#^^' 



b^a 



b^a 



+ h.ylya [7eAaaoM<, + O-ao^aaoX + 2aaoe«aoA + 2Wao£«aoA + 7eA E ^^ + ^ E 



diubx . dwbe 



b^a 



b^a 



dy'a dy^a 



+ ro«eyl)iya^yil)} + o{c"% 
(ya) - ~7o.x^{Y.^^(y-)-T.{'^b+yl^)'To{^yayl~i''ya,y!i}<e} + o{c-'), 

^ b b^a ^" 

9a(3iya) = laP+lal3-^[^Wb{ya)-^[Wb+yl-^jj+0{c-'^), 



90 



b^a 



with the equation of motion of the body a, x^o = ^"o + ^ ^Qao + ^(^ ^) given by Eq. ([236]). 

The coordinate transformations that place an observer into this reference frame are given below: 



„o _ „,o 



ya+C '{y"^°(E"'^~lK.+C ^qaoe)K^+C~^qlj)dy'^'-civ^^,+C ^aoe)yl} 

yan b^n 



"0 b^a 



+ c- 



l(y y^ )2_ 3 



7Ap 



^K^ , d^bx 



b^a b^a b^a 



aoe^ao "•" '^aoX^aofi. "•" '^aofi^aoX ' 



E^^ + 25:(: 



b^a 



b^a 



dya dya 



a + ^CpM^', 



(251) 
(252) 



(253) 
(254) 

(255) 



(256) 



^" - ya+<+c~' 



{C - (^ao<„ + 7"^ E ^^ + <) y'^^ + «aoe (y" 2/: - h^'y^A^a) } + 0{c-^). (257) 



b^^a 



The expressions (j253|) - (|255p represent the harmonic metric tensor in the local coordinates of the accelerating 
reference frame. The form of the metric (|253p - (|255p and the coordinate transformations (|256p - (|257p are new and 
extend previous formulations (discussed, for instance, in [ll|) obtained with different methods. These results can 
be used to develop models for high precision navigation and gravitational experiments. However, for a complete 
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description of these experiments we would need to establish the inverse coordinate transformations - the task that 
will be performed in the next section. 

Finally, we note that the results presented here can be generalized to extended bodies. After a tedious but straight- 
forward calculation (which can be done in a manner similar to [25|), we find that, given the accuracy of the measure- 
ments in the solar system, only the following two modifications need to be made to our solution: 

ICaiVa) - l\Y.(^''^a~h^aoyao)dya-civ,^eyl) + Oic-^), (258) 

•'yao b^a 

< = -Y.(9''Ub), + Oic-*). (259) 

The difference between expressions for JCaiUa) given by Eqs. (|220p and (|258p and also between those derived for a" 
and given by Eqs. (|22ip and (|259p is in the fact that the new expressions account for the interaction of the body's 
gravitational multipole moments with the background gravitational field, via the procedure of averaging the external 
gravitational field over the body's volume in the form of 

{Ub)^ = m-i / d?y',aUb, {d^Ub}, = m^' f d'y'.ad^Ub, m^ = f d^y'^a, (260) 

as opposed to accounting only for the value of that field on the body's world-line (via the limiting procedure defined 
by Eq. (PTU)) taken for Ub and d"Ub). 

V. INVERSE TRANSFORMATIONS 

In the preceding sections, we constructed an explicit form of the direct transformation between a global inertial and 
local accelerating reference frames by applying the harmonic gauge and dynamical conditions on the metric. Given 
the Jacobian matrix Eqs. (|42p - (|43p . it was most convenient to work with the covariant form of the metric tensor, 
which could be expressed in terms of the accelerating coordinates by trivial application of the tensorial transformation 
rules Eq. ([T^ . 

The same logic suggests that if we were to work on the inverse transformation: that is, when it is the inverse 
Jacobian matrix {dy"^ / dx"} that is given in explicit form, it is more convenient to work with the contravariant form 
of the metric tensor, to which this Jacobian can be applied readily. This simple observation leads us to the idea 
that we can get the inverse transformations — i.e., from the accelerated to the inertial frame — by simply repeating the 
previous calculations, but with the contravariant form of the metric tensor instead of the covariant form. 

In this section, we show that this is indeed feasible, and accomplish something not usually found in the litera- 
ture: construction of a method that can be applied for both direct and inverse transformations between inertial and 
accelerating reference frames at the same time, in a self-consistent manner. 

A. The contravariant metric tensor in the local frame 

We have obtained the solution of the gravitational field equations for the A^-body problem, which in the barycentric 
reference frame has the form Eqs. (|116p - (|118p and (|119p . (|120p . It follows from Eq. ([5T|) that the solution in the local 
frame could be in the following form: 

9-{Va) = ^^9'\<ya)). (261) 

Using the coordinate transformations given by Eqs. (P(7)) - (|^T|) together with Eqs. P9|) - ([5(I)) . we determine the form of 
this metric tensor in the local reference frame associated with the body a: 

00/ ^ 1 , '2' (dta , 1 ^x'^dtaldtaX , 

^ ^y^^ = i + ^la^ + 5^ -c^-c^\^ 

c* I dx'^ ^ V dx'^ / c dx'^ c dx^ V dx'^ ^ c dx'^ c dx^ J J 
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+ ^(E-^(-(y»))+a^ + ^7"-a;r-a^) +'^(c^^)^ (262) 

b 

+ ?r c9^ + '^9^ + ^ ca^a^-"-9^l + '^(^ )' (263) 

5"/5(y,) = ^'^P + ^{y2yao+l"'^+^''^}-r'^j:^>.Ky.)) + Oic-% (264) 

b 

As in the case of the direct transformation, we impose the harmonic gauge condition on the metric, to help us 
estabUsh cxpHcit forms of the transformation functions JCa, Ca, and Q". 

Using Eqs. ((80l) and ((82)) to constrain the form of the metric tensor, given by Eqs. (|262p ~ (|264p in the local frame, 
we obtain: 



2 NT- r/. 2 , \ , , ,, 4 

b 

r{ya) = |E(<(-(2/'^))-<o-^(-(2/a)))+^{7"^^^ + c^+<og-<|^}+0(c-^), (266) 



b 



^iVa) = 7"^ + ^{<o<+7"^|^+7^^^}-7"^|E-''(-(J/'^))+^(^-')- (267) 

b 

In analogy with Eq. (|96l). we determine the source term, S"^^{ya), for the A^-body problem in the local reference 
frame as: 

S""\ya)-^^S^\xiya)), (268) 

where 5^'(a;(j/a) is the source term in the global frame. Thus, S^"'"{ya) takes the form: 

S°'{ya) ^\c^ E { (l - |«ao.<c)^^(^(ya)) + ^^^ao.^^^ (^(^a)) + |(|^ + ^ao.^o) '^"(^(y-)) + ^(c"')}^ (269) 

^°"(ya) =c^ {a?(a:(2/a)) - t':„ab(x(2/a)) + 0(c-^)}, (270) 

^"^(ya) = - y\^ E {'^^(^(ya)) + 0{C-^)]. (271) 

6 

The metric tensor ©-(J!]) and the partial-differential form of the gauge conditions ([73|) -([74 | together with condition 
(15^ allow one to simplify the expressions for the Ricci tensor and to present its contravariant components i?™" — 
g"^''g"^Rki in the following form: 

R'"{ya) = in,„(c-2.g[2l0° + c-^{.gW00-l(g[2l0")^})+O(c-6), (272) 

R^'^iVa) = c-3iA,„5pl0" + O(c-5), (273) 

R'^^'iVa) = c-21a,„5[21"^ + 0(c-4). (274) 

These expressions allow one to express the gravitational field equations as below: 

2 , \ , . .. 4 



b 
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A-kG 



2 .a^a , 1 



E{(l-^"aoXo)'^''(^(2/a)) + ^«ao.'^^(^(y«)) + ^(5^ + l«ao.<o)^^(^(ya))+0(c^')}, (275) 



b 

= 4^G5]{a?(x(yJ) -<„ab(x(yj)} +0(c-2), 



(276) 



1 / (90^ BO°'\ n 

y"'^5^u;,(x(y,))--(««^f„ + ^"^^+/^^)+0(c-2)J =4^G7"^^afc(:c(y„))+0(c-2)(.277) 



Remembering Eq. (|150p and noting that lS.y^ = A-c + ©(c ^), with help of Eqs. ([86]) - ((88)) . we can verify that 
Eqs. (P75)) - ((7r7)) are satisfied. 



m /^„,fc 



We now note that the Jacobian matrix on the right hand side on Eq. (j26ip is composed of the quantities dy"' /d. 
that arc functions of {x''}, the same is true for g*'^{x{ya)): which may also be treated only as a function of {2:*''}. 
Therefore, the metric on the right hand side can also be treated as a function of g™-"{ya{x)): 






(278) 



with the understanding that (/'""(a;) — g'""''{ya{x)). Components oi g"''^{x) can be expressed, similar to Eqs. (|15ip - 
(jl53p . using the two harmonic potentials: 



gOO(x) = l + \njix) + ^W^ix)+Oic-'), 



g^><-{x) = u;"(x)+0(c-^), 

g^P(x) - 7"^_7"/^l^(x-)+0(c-^), 



(279) 
(280) 
(281) 



where the total (gravitation plus inertia) local scalar 'w{x) and vector 'w°'{x) potentials have the following form: 



w{x) = y^ Wfc(x) 



dx° 



+ ^v„^,v„^ + 



1 (dCa , vl^^dCa , i/a^a\2 /S/C^ , ^ 



2 aoe ao ' „2 



la^ + Va^ + Ha^) -(a^+2-aoe-aoj }+0(c ),(282) 



W 



'{x) = ^<(x) + i{ 






}+0(c-'), 



(283) 



with gravitational potentials in the local frame Wb{x) and w^lx) expressed in global coordinates related to their 
counterparts in the global frame 'Wb{x{ya)) and w^{x{ya)) as below 



M^) = (1 - -^v^^y^^jwtix{ya)) + -^v^^^wlixiya)) + 0{c '^), 
w^{x) = wnx{ya))~vZ,w,{x{ya))+0{c~^). 



(284) 
(285) 



Clearly, the potentials Wb{x) and w^{x) introduced by Eqs. (j282l) - (|283p satisfy the following harmonic equations 
in the local coordinates {j/™}: 



2 ,dK„ 



UyMx) = 4^G^ {(1 - -v,^^,vl^yb{x{ya)) + ^v,^,al{x{ya)) + ^(^ + \va,eVl,)^h{x{ya)) + 0(c-4)|, 



\yw"{x)=A^GY,[a'S{x{ya))-v:M^iya))+0{c-^)]. 



(286) 
(287) 
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With the help of Eqs. (I86l)-(l88l) one can verify that both Eqs. ((286)) and ((287)) are identicahy satisfied. 

Finally, using the relations cd/dy^ = cd/dx° + vl^d/dx^ + ©(c'^) and d/dy^ = d/dx"" + ©(c^^) in Eq. (gl]), one 
can verify that potentials Wb{x) and w^{x) satisfy the following continuity equation: 

d . d \ , , d 



These expressions allows us to present the scalar and vector potentials Eqs. ((328)) ~ ([329]), expressed entirely as 
function of the global coordinates. Remembering Eq. (|147p and also the fact that 

a^Wb{x{ya)) ^ 47rGabix{ya)) + 0{c-^), A^w^ix{ya)) = A7rG<j^{x{ya)) + 0{c-^), (289) 

where solutions are given as: 

w,{xiya)) = G /dV^|iM^+ 1 G^ /dVa,(2/,(x'))|x-x'|+0(c-4), (290) 

J |x-x'| 2c^ dx^ J 

w^xiVa)) = G jd^x'^^^+0{c-% (291) 

we can use these results to present the gravitational potentials (|290P " ((29ip in coordinates of the global frame in the 
form of the integrals over the body's volume as below (where x™ is understood as a;™ = x™'{ya)): 

+ ^G^y'dVa,(2/,(x'))|x-x'|+0(c-3), (292) 

,S.i^b{yb{x')) „,a r^ [ .3,i^b{yb{x')) , ,^. -2n 



<(x) =GJ d'x'^^^^ - vlfij d'x'^^^^ + 0{c~% (293) 

Although the expressions for the scalar inertial potentials have different functional dependence on the transformation 
functions (i.e., {IC,C, Q") vs. {IC,C, Q")), it is clear that both expressions for w{ya) given by Eqs. (|154p and ((282|) 
are identical. The same is true for the inertial vector potentials w°'{ya) given by (|155p and (|283p . 

1. Determining the structure of Ka 
The general solution to Eq. ((55)) can be thought in the following form: 
Ka{x) = ka,+ka,^ri^a+6ka{x)+0{c-^), whcrc '5Ka(a:) - ^ -^ Ac,„^, ...^, (rO)r^i -T^'^ + ^(c-^) , (294) 

where r" is defined as below: 



A:! 

fc=2 



C = -^^ - <o . <, = < + c^'C, + 0{cr% (295) 

where a;"^(a;*') is the position vector of the body a, complete to 0{cr^), and expressed as a function of the global 
time-like coordinate x^ , given by Eq. ((47)) and ^^^^^...^^(a;*') being Cartesian spatial trace-free (STF) tensors [23| 
which depend only on the time taken at the origin of the coordinate system on the observer's world-line. Substituting 
this form of the function ICa in the equation ((9T)). we find solution for kaa^ and Ska'- 

kaot. = ™aoM + 0(C"^), 5ka = ©(c'^). (296) 

As a result, the function Ka that satisfies the harmonic gauge conditions was determined to be 

ta{x) = kao + c(Vao,.r^) + 0{c-^). (297) 

This expression completely fixes the spatial dependence of the function ICa , but still has unknown dependence on time 
via function kaa{x'^). This dependence determines the proper time on the observers world- line which will be done 
later. 
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2. Determining the structure of Q" 

The general solution for the function Q^ that satisfies Eq. ((88)) may be presented as a sum of a solution of the 
inhomogcncous Poisson equation and a solution of the homogeneous Laplace equation. Furthermore, the part of that 
solution with regular behavior in the vicinity of the world-line may be given in the following form: 

Q:i^) = -< + €„,< + ^€oM.« + C(^) + 0{c-'), (298) 

where g" „j, can be determined directly from Eq. (|85p and function S^" satisfies the Laplace equation 

We can see that Eq. (|88p can be used to determine qaain^i but would leave the other terms in the equation unspecified. 
To determine these terms, we use Eq. (|92|) (which is equivalent to Eq. ([88]) ) together with Eq. (|9T|) . and obtain: 

«ao«ao+7 Q.^X+1 dx^ ^ V 9^0 ^ "'^o.^aoj ^l^ ). (6VV ) 

Using the intermediate solution for the function fCa from Eq. (|297p in Eq. (|300p . we obtain the following equation 
for Q«: 

<-fo+7"^gf + /^|^ + 27"^(|:^-KoXo+«oe<) = Oic~-). (301) 

A trial solution to Ec^. pOip may be given in the following general from: 

-u;:X,+6t{x)+0{c--'), (302) 

where q" and antisymmetric matrix w"^^ = — w^" are some functions of time, and ci,...,C5 are constants and 5^l^{x) 

is given by Eq. (|299p and is at least of the third order in spatial coordinates r^, namely S^l^{x) oc Odr^'j^). Direct 
substitution of Eq. p02p into Eq. poip results in the following unique solution for these coefficients: 

d = --, C2 = 0, C3 = -, C4 = -1, C5 = -1. (303) 

As a result, function Q" has the following structure 

Q:i^) = -C„-(Ko<+<+7"^(|^-HoA<))rae-«.„e(^X-h"XA^a)+C(^)+e^(c"'), (304) 

where g"o ^'^^ '^"o ^^^ y^^ ^'^ ^e determined. 

By substituting Eq. p04p into Eq. pOip . we see that the function S^2{x) in Eq. p04p must satisfy the equation: 

A-jS^^^SC^Olc-'y (305) 

We keep in mind that the function 5^^ (x) must also satisfy Eq. (|299p . The solution to the partial differential equation 
(|299p with regular behavior on the world-line (i.e., when |ra| — > 0) can be given in powers of r^' as 

^iaix) =Y.h^^ao,r...,A^°K'-<' +0{K\^)+0{C-'), (306) 

fc>3 

where S£,2 uj Ufc(^°) being STF tensors that depend only on time- like coordinate xq. Using the solution p06p in 
Eq. pOSp . we can see that S^^ u u is a-lso antisymmetric with respect to the index a and any of the spatial indices 
/ii.../ifc. Combination of these two conditions suggests that <5Caoui u„ ~ ^^r all fc > 3, thus 

di:{x)=0{c-'). (307) 

Therefore, application of the harmonic gauge conditions leads to the following structure of function Q": 

Q^ai^) = -Co - (l<o< +< +7"1^ - HoxOYae ' %,. (C< " b'^'r^.V^) + 0{c''), (308) 

where Coi"^™ ^'^'^ '^"■a '<^^^ yet to be determined. 
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3. Determining the structure of La 



We now turn our attention to the second gauge condition on the temporal coordinate transformation, Eq. (|87l 
Using the intermediate solution (|297|) for the function JCa, we obtain the following equation for Ca'- 






The general solution of Eq. p09p for £a may be presented as a sum of a solution 5£a for the inhomogeneous Poisson 
equation and a solution 6CaQ of the homogeneous Laplace equation. A trial solution of the inhomogeneous equation 
to this equation, SCa, is sought in the following form: 

SU^) = cfciK,<J(r,^<) + cfe«,r:)(r,,<) - fcgc^ A(^ + i«,^,<,) (r,,<) + Oic-% (310) 

where ki,k2,k3 are some constants. Direct substitution of Eq. piOp into Eq. p09p yields the following values for 
these coefficients: 

111 

ki = g, fca = — ' ^3 = g- (311) 

As a result, the solution for 6Ca that satisfies the harmonic gauge conditions has the following form: 



6C. 



a{x) = ic(^;,^,<J(r,,<)-ji^c«,r:)(r„,<)-ic2^(^ + i^;,„,<J(r,X) + 0(c-')• (312) 

The solution for the homogeneous equation with regular behavior on the world-line (i.e., when |r^| — > 0) may be 
presented as follows: 

>Cao(x) = £aoix°)+Lox{x'')r^ + ^£aMx'')r^r^ + Siaix) + Oic-^), (313) 

where iao\fi is the STF tensor of second rank and 6ia is a function formed from the similar STF tensors of higher 
order 

SL{x) = Y.^_SL„,,...,,{x'K^-rii''+0{Kn+0{c-'). (314) 

fe>3 

Finally, the general solution of the Eq. (|309p may be presented as a sum of the special solution of inhomogeneous 
equation SCa and solution, Cao, to the homogeneous equation, ACa = 0. Therefore, the general solution for the gauge 
equations for function Ca (y) = Cao + 6Ca has the following form: 

Calx) = lag + «aoA ''a ' 2 '^oAp^a'^a + ^'^y^aot'^ao) Vai^'^'a) ~ To'^\'^aoe^a)Vai^'''a) ~ 

- I^"^{^ - h'aoXo)ira.r:)+SL{x)+0{c-'). (315) 

We successfully determined the structure of the transformation functions /Ca, Q" and Ca, which are imposed by 
the harmonic gauge conditions. Specifically, the harmonic structure for K. is given by Eq. (|296p . the function Q" was 
determined to have the structure given by Eq. p09p , and the structure for Ca is given by Eq. pisp . The structure of 
the functions JCa, Q" and Ca given in the following form: 

ICaix) = ka„+ciVa„^r^)+0{c-''), (316) 

Q:{^) = -Co - {i<Ao + < + ^"'(^ - baoxVa„))rae - aa„,{ryi - h^'raxr',) + 0{c-'), (317) 

Ca{x) = iao + iaoX r^a + ¥a„X^. r^r^ + ic(2(z;,„,<J - ''d^^^ " h'aoeO) ('^aA^a) " 

- ^cia,^yaKraXa) + 6iaix)+Oic-'). (318) 

Note that the harmonic gauge conditions allow us to reconstruct the structure of the functions with respect to spatial 
coordinate x'*. The time-dependent functions k^q , g"^ , w"^*^ , ^aoj^aoA,^aoA/ii and (^^ao^i...^^ cannot be determined from 
the gauge conditions alone. Instead, just as we did in Sec. IIVI we resort to a set of dynamic conditions to define 
unambiguously the proper reference frame of a moving observer. 
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B. Finding the form of the coordinate transformation functions 

111 the case of the metric tensor given by the expressions Eqs. (|279p - (|28ip . the conditions Eqs. (|208p - (|209p with 
r^ = x"^ — x^^j given by Eq. (|295p lead to the following set of equations. 

lim wix) = Wa{x)+Oic-% lim ^^ = ^^ + 0(c-4), (319) 

|r„|-s-0 |r„|-i.O OXP OxP 

lim u;"(x) = <(x) +0(c-2), lim ^^ = ^^+0{c~'). (320) 

Ir^l^O |r„|-i.O OxP OXP 



Imposing the conditions given by Eqs. pi9|) - p20p on the potentials w and w" , which are in turn given by Eqs. ()282p - 
(|283p . results in the following set of partial differential equations set on the world- line of the local observer expressed 
via global coordinates {x™}: 






^^b + -Q^ - ^V^oeVao ' ^ "Wao.Co 






Z.5^+°ao/J + ^|^^^^ + — ^^^^-aao/3^^1 - ^l^ j, (622) 

"'^+n^ c9^ + ^9^ + --a^-"-(a^---0| = 0(c ), (323) 



with all the quantifies above now being the functions of the global time- like coordinate x'^ . 

The equations above can now be used to determine uniquely the form of the coordinate transformation functions. 
From the first two conditions p2ip and (|322p above, we immediately obtain 

E t^^ + ^ - ^-aoe<„ - c-\gj:, = 0{c-% (325) 

fill 

Efe|+«ao^ = 0{c-% (326) 

where we split the scalar potential Wb{x) defined in Eq. (|284p into Newtonian (Ub) and post-Newtonian {6wb) parts: 

Wb = Ub + \5wb + 0{c-^), (327) 

so, that gravitational scalar and vector potentials defined by Eqs. (|284p - (|285p take the following form: 

Wb{3)a{x)) = Ub{x) + 1 {Sw,{x) - 2{Vg^^yjUb{x) + 4v,^wt{x)) + 0{c-^), (328) 

wUVaix)) = w-S{x)-vlPb{x)+0{c-^). (329) 

Eqs. (|325p and (|297p allow us to determine function kg^ and to present the solution for K-a as below: 

^a{x) = - [Y.Ub-\v^,,vl^-c-\^Jl^)dx'^ + c{v,^jl) + 0{c-^). (330) 

■^^S b^a 

Substituting this expression into Eq. pi7p . we can determine Q": 

QZ{X) = -Co + (7"^ Y.^b- i<,<o - <)''- - aaoe(C< " h'^'^aXT^) + ^(c"')- (331) 

b=^a 
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Finally, the general solution, given by Eq. (|315|) . for the function Ca takes the following form: 



Ca{x) = L, + iaoA r^ + \ (^oA,. + C7A,. (f ("aoeO + \c-^ ^ t/fc)r^< - l^c(a,„,<)( 



auaJ ' 



-2\ 



+ ^4(x)+0(c-"). (332) 

This equation allows us to proceed onto finding the function Ca- To do this, we rely on the yet-unused parts of 
Eqs. (|32T|) -((324 l) which take the following form 

E^-^+la^ + V^^ + Hs^-^-^-J - ( a^ - 2-aoe-ao) | = ^(c ), (333) 



^< + H^"^^|^ + ^^ + <S|-<(^-W:j} - Oic-% (335) 



ia£ 



Using the expressions for the functions K-a, Q" given by Eqs. p30p and p3ip correspondingly together with the 
expression p32p for £a, and also using transformation rule for the gravitational potentials given by Eqs. p28p and 
P29p together with the Newtonian equations of motion p26p . we can now obtain the equations needed to determine 
the remaining unknown functions £ao,4oA, ^aoAp, ^ao ^'^'^ ^Zr 

Eq. p33p leads to the following equation for ia^,- 

Next, Eq. f|334p results in the equation for l'^^: 

-Xo = - E ^J^" - U^ao.<oX - <oE ^b - 4^ao. E ^"'^^ + ^(^"')- (338) 

b^a b^a b^a 

From Eq. (|335p we determine i^^ : 

-Jao - Co - Ho i^aoeVl ) + 3^ E ^^ " ^ E ^'" + ^(c"')' (339) 

Furthermore, Eq. (|336p leads us to the following solution for £"^^: 

b^a ^ ^T^a 

Using the same argument concerning the symmetry properties of i^o that led to Eq. p32p , we find the following 
unique choice for the anti-symmetric matrix Coa,, that symmetrize Eq. p40p : 

b^a ^ 

which, if not accounting for Newtonian equation of motion f|326p . actually has the following form: 

^:f = -U<<-<yaJ-^T.{<§^-<^)-^j:i9-<-9'<)+0{c-^). (342) 
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The first term in the equation above is the special relativistie Thomas precession, the second term is the geodetic 
precession and the last term is the Lense-Thirring precession. 

The result p4ip allows us to present Eq. p40p for the function ^"^ in the following form: 

^Cf = -2E(U + S)-|7"^E#-K«fo+-fo< + |7"^..Xo)+e)(c-). (343) 

Finally, Eqs. p38p and ([339]) allow us to determine the equation for q^^. Indeed, differentiating Eq. (|339p with 
respect to time and subtracting the result from Eq. (j338p . we obtain: 

h^a b^a b^a b^a 

f)Tj 
-3<o E ^ai - (^'.o.^oXo + <„(-ao.<o) + 0{c~^). (344) 

By combining Eqs. p26p and (|344p . we obtain the equations of motion of the body a with respect to the inertial 
reference frame {x™}: 

b^a b^a b^a b^a 

dtJb 
The equations of motion p45p establish the correspondence barycentric equations of motion of the body a. 



3<oE^|^ - K.oe<o)<o +CKe<o)) +e^(c-^)- (345) 



C. Summary of results for the inverse transformation 

In this section, we sought to write the transformation between global coordinates {x^} of the barycentric reference 
system and local coordinates {y^} introduced in the proper body-centric coordinate reference system associated with 
the body a in the form: 

yl = X° + C^',^a(x'=)+C-4£a(2:'=)+0(c-^), (346) 

y: - x" - zl^{x^) + c-^QZix"") + 0(c-4). (347) 

Imposing the same conditions discussed in Sec. llVCi we were able to find explicit forms for the transformation 
functions Ka , C,a and Q" : 



^a(x) = -l[Y.U,- ^^v,^Xo - c-'vaoMo)d^" + c{v^,A) + 0{c~^), (348) 

Sa(^) = -ra„ + {r'Y.Ub~\<<o-<:y<^^-''<^^e{<rl-\r'ra>,rl)+0{c-^). (349) 

were the anti-symmetric relativistie precession matrix uj'^^ in the barycentric reference frame is given as below: 

<^ = li« - «) - 2 E (U - S) + ^(^-^)- (350) 

b^a "^ 

Also, the function Ca was determined in the fohowing form: 

Ux) = 4„+^;„,r,^-icr,,r,,[<a(^„-f<,a,^„+7'^E^^ + 2E(|^ + |^ 

b^a b^a ^ 



38 
-To«er:)ira.r:) + 5L{x) + 0{c-% (351) 

with functions lao ^-nd f^^ are given by 



b^a 



b^a 



b^a 



b^a 



-Ao = Co-KoKe<o) + 3<E^''-4E<+^(c"')- 



(353) 



6#a 



b^a 



Finally, the function q"^^ is the post-Newtonian part of the position vector of the body a given by Eq. (|344p . Together 
with the Newtonian part of the acceleration, Za„ = a" , the post-Newtonian part q^^ completes the barycentric equation 
of motions of the body a given by Eq. ([345]). This equation essentially is the geodetic equation of motion written for 
the body a in the global reference frame. 

Substituting these solutions for the functions ICa , Q" and Ca into the expressions for the potentials w and w" given 
by Eqs. (|282p - (|283p . we find the following form for these potentials: 



\\2 a' a 



b^a 
1 
72 



3aao£«aoA + «Qoe«QoA + '"a„eO'ao\ + 



T.aE^^ + ^E 



b^a 



b^a 



dx' 



^6A 



dx-^ 



•^bt 



w 



"(-) = T.<(-^ -T.{< +<^) - To{^<<-r'vii}<. + Tc^^^a+o{c- 



(354) 
(355) 



b^a 



The same argument that allowed us to eliminate 5£a in Sec. IIV CI works here, allowing us to omit 5ia- Therefore, 
we can now present the metric of the moving observer expressed in the global coordinates {x"}. Substituting the 
expressions for the potentials given by p54p and p55p into Eqs. p79p ~ (|28ip we obtain the following expressions 



b^a 



b^a 



i<^a(3aao.«aoA+a 



aoe aoA "*" ^ a^e a^X 



+ l^\^Ub + 2 ^(a^WbA + ^AWfc) 



U<oya){ra,<)]+O{c-^). 



b^a 



b^a 



dw? 



5°"(^) = ;3{E<(^)-E(<+<^)-ii){3C'^^-T"XX}<^}+^(^'')' 






b b^a 

2 



{^.,(.)-5:(^,+.:||)}+o(c-). 



(356) 
(357) 

(358) 



b b=^a 

The coordinate transformations that put the observer in this reference frame are given by: 



Vl = -° + c-' 



+ C [{Vaoe + C-Haoe) (l + C'^ (3 E ^^ " l^aoA) ^ c"'4^ U;,,) r^ } + 
•'^o b^a b^a b^a 

\cr,xr,^(yla':,^,+vlal+^>^^Y.^, + 2Y,{d^w>i+d^wl)) - i^c«,r:)(r,X)} + 0(c-'), (359) 



65^0 



b^a 
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The obtained inverse coordinate transformations are new and extend the previously obtained results (for review, 
see [5|, 0], and also [gl). This set of results concludes our derivation of the coordinate transformations from a local 
accelerated reference from to the global inertial frame. Generalization to the case of N extended bodies can be done 
in a complete analogy with Eqs. ([258]) -([260]), which is sufficient for the solar system applications. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we introduced a new approach to find a solution for the A^-body problem in the general theory 
of relativity in the weak-field and slow motion approximation. The approach is based on solving all three problems 
discussed in Sec. HI including finding solutions to global and local problems, as well as developing a theory of relativistic 
reference frames. Our objective was to establish the properties of a local frame associated with the world-line of one 
of the bodies in the system and to find a suitable form of the metric tensor corresponding to this frame; wc achieved 
this goal by imposing a set of clearly defined coordinate and physical conditions. 

Specifically, we combined a new perturbation theory ansatz introduced in Sec. Ill CI for the gravitational A^-body 
system. In this approach, the solution to the gravitational field equations in any reference frame is presented as a sum 
of three terms: i) the inertial flat spacetime in that frame, ii) unperturbed solutions for each body in the system boosted 
to the coordinates of this frame, and iii) the gravitational interaction term. We use the harmonic gauge conditions that 
allow reconstruction of a significant part of the structure of the post-Galilean coordinate transformation functions 
relating non-rotating global coordinates in the inertial reference frame to the local coordinates of the non-inertial 
frame associated with a particular body. The remaining parts of these functions are determined from dynamical 
conditions, obtained by constructing the relativistic proper reference frame associated with a particular body. In this 
frame, the effect of external forces acting on the body is balanced by the fictitious frame-reaction force that is needed 
to keep the body at rest with respect to the frame, conserving its relativistic linear momentum. We find that this 
is sufficient to determine explicitly all the terms of the coordinate transformation. Exactly the same method is then 
used to develop the inverse transformations. 

Our approach naturally incorporates properties of dynamical coordinate reference systems into the hierarchy of 
relativistic reference systems in the solar system and relevant time scales, accurate to the c~^ order. The results 
obtained enable us to address the needs of practical astronomy and allow us to develop adequate models for high- 
precision experiments. In our approach, an astronomical reference system is fully defined by a metric tensor given 
in a particular coordinate frame, a set of gravitational potentials describing the gravitational field in that frame, a 
set of direct, inverse, and mutual coordinate transformations between the frames involved, relevant time arguments, 
ephemerides, and the standard physical constants and algorithms. We anticipate that the results presented here may 
find immediate use in many areas of modern geodesy, astronomy, and astrophysics. 

The new results reported in this paper agree well with those previously obtained by other researchers. In fact, the 
coordinate transformations that we derived here are in agreement with the results established for both direct and 
inverse coordinate transformations given in Refs jlll - [l4l | and [a, [S Ej [l3 1 correspondingly. However, our approach 
allows one to consistently and within the same framework develop both direct and inverse transformations, the 
corresponding equations of motion, and explicit forms of the metric tensors and gravitational potentials in the various 
reference frames involved. The difficulty of this task was mentioned in Ref. [y] when the post- Newtonian motion of a 
gravitational A^-body system was considered; our proposed formulation successfully resolves this important issue. As 
an added benefit, the new approach provides one with a good justification to eliminate the functions 5k, 5^ and 51, 
yielding a complete form for the transformation functions /C, C and Q" involved in the transformations (as well as 
their "hatted" inverse counterparts). 

The significance of our result is that for the first time, a formalism for the coordinate transformation between 
relativistic reference frames is provided, presenting both the direct and inverse transformations in explicit form. By 
combining inverse and direct transformations, the transformation rules between arbitrary accelerating frames can be 
obtained. Furthermore, it is possible to combine direct (or inverse) transformations, and obtain a complete set of 
transformations that can be represented by our formalism, as shown explicitly in [8|, ll6| . This leads to an approximate 
finite group structure that extends the Poincarc group of global transformations to accelerating reference frames. 

The results obtained in this paper are designed to facilitate the analysis of relativistic phenomena with ever increas- 
ing greater accuracy. We should note that the approach we presented can be further developed in an iterative manner: 
if greater accuracy is desired, the coordinate transformations ([57)) - ((55|) can be expanded to include higher-order terms. 
Furthermore, the same approach relying on the functional /Ci2Q-parameterization may be successfully applied to the 
case of describing the gravitational dynamics of an astronomical A'^-body system and dynamically rotating reference 
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frames. This work has begun and the resuhs, when available, will be reported elsewhere. 

There are some problems that remain to be solved. One of these is the problem of relativistic rotation. In particular, 
it is known that the rotational motion of extended bodies in general relativity is a complicated problem that has no 
complete solution up to now. On the other hand, modern observational accuracy of the geodynamical observations 
makes it necessary to have a rigorous relativistic model of the Earth's rotation. The method presented in this paper 
may now be used to study the rotation of extended bodies to derive a theory of relativistic rotation within the general 
theory of relativity. Our ICCQ formalism presented here provides one with the necessary conceptual basis to study 
this problem from a very general position and could serve as the foundation for future work. Also, a complete formal 
treatment of the case of a system of N extended bodies should be presented. Although this treatment would result in 
corrections much smaller to those anticipated in the solar system experiments, they may be important to a stronger 
gravitational regime, such as in binary pulsars. This work had been initiated. The results will be reported elsewhere. 

Our current efforts are directed towards the practical application of the results obtained in this paper (similar 
to those discussed in [26| in the context of the GRAIL mission). We aim to establish the necessary relativistic 
measurements models, equations of motion for the celestial bodies and spacecraft, as well as the light propagation 
equations in different frames involved; we also plan to implement these results in the form of computer codes for the 
purpose of high-precision spacecraft navigation and will perform the relevant scientific data analysis. The analysis of 
the above-mentioned problems from the standpoint of the new theory of relativistic astronomical reference systems 
will be the subject of specific studies and future publications. 
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